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ABSTRACT 

Using 7 years of MACHO survey data, we present a new determination of the optical depth to mi- 
crolensing towards the Galactic bulge. We select the sample of 62 microlensing events (60 unique) on 
clump giant sources and perform a detailed efficiency analysis. We use only the clump giant sources 
because these are bright bulge stars and are not as strongly affected by blending as other events. Using a 
subsample of 42 clump events concentrated in an area of 4.5 deg 2 with 739000 clump giant stars, we find 
r = 2.17to.tl x 10_6 at (') b ) = ( 1?50 ' -2?68), somewhat smaUer than found in most previous MACHO 
studies, but in excellent agreement with recent theoretical predictions. We also present the optical depth in 
each of the 19 fields in which we detected events, and find limits on optical depth for fields with no events. 
The errors in optical depth in individual fields are dominated by Poisson noise. We measure optical depth 
gradients of (1.06 ± 0.71) x lO^deg -1 and (0.29 ± 0.43) x lO^dcg -1 in the galactic latitude b and 
longitude I directions, respectively. Finally, we discuss the possibility of anomalous duration distribution 
of events in the field 104 centered on (l,b) = (3? 11, — 3?01) as well as investigate spatial clustering of 
events in all fields. 

Subject headings: Galaxy: center, stellar content, structure — gravitational lensing — surveys 
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1. Introduction 

The structure and composition of our Galaxy is 
one of the outstanding problems in contemporary as- 
trophysics. Microlensing is a powerful tool to learn 
about massive objects in the Galaxy. The amount 
of matter between the source and observer is typ- 
ically described in terms of the microlensing opti- 
cal depth, which is defined as the probability that a 
source flux will be gravitationally magnified by more 
than a factor of 1.34. Early predictions (Griest et al. 
1991; Paczy hski 1 99 1 ) of the optical depth towards the 
Galactic center included only disk lenses and found 
values near t = 0.5 x 10~ 6 . The early detection 
rate (Udalski et al. 1993, 1994a) seemed higher, and 
further calculations (Kiraga & Paczyhski 1994) added 
bulge stars to bring the prediction up to 0.85 x 10~ 6 . 
The first measurements were substantially higher than 
this: t > 3.3 ± 1.2 x lO" 6 at (l,b) w (l?0,-3?9) 
(Udalski et al. 1994b) based upon 9 events and r = 
3.9t{- 8 2 x 10~ 6 at (l,b) = (2?52,-3?64) (Alcock 
et al. 1997a) based upon 13 clump-giant events and 
an efficiency calculation. Many additional theoret- 
ical studies ensued, including additional effects, es- 
pecially non-axisymmetric components such as a bar 
(e.g., Zhao, Spergel & Rich 1995; Metcalf 1995; Zhao 
& Mao 1996; Bissantz et al. 1997; Gyuk 1999; Nair 
& Miralda-Escude 1999; Binney, Bissantz & Ger- 
hard 2000; Sevenster & Kalnajs 2001; Evans & Be- 
lokurov 2002; Han & Gould 2003). Values in the range 
0.8 x 10~ 6 to 2 x 10~ 6 were predicted for various 
models, and values as large 4 x 10~ 6 were found to be 
inconsistent with almost any model. 

More recent measurements using efficiency calcu- 
lations found values of r = 2.43l°;38 x 10 ~ 6 at 
(I, b) = (2?68, -3?35) from 99 events in 8 fields using 
difference image analysis (Alcock et al. 2000a; MA- 
CHO), t = 2.59toit x 10~ 6 at (I, b) w (1?0, -3?9) 
from 28 events using difference image analysis (Sumi 
et al. 2003; MO A) 18 , r = 2.0 ± 0.4 x 10~ 6 at 
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We would like to notice here that the results from the difference 
image analyses cannot be easily compared to the ones from clump 
analyses. The analyses that use clump giant sources have reasonably 



(l,b) = (3?9,-3?8) from around 50 clump-giant 
events in a preliminary version of this paper (Popowski 
et al. 2001a; MACHO), and r = 0.94 ± 0.29 x 10~ 6 
at (l,b) = (2?5,-4?0) from 16 clump-giant events 
(Afonso et al. 2003; EROS). 

Blending is a major problem in any analysis of the 
microlensing data involving point spread function pho- 
tometry. The bulge fields are crowded, so that the 
objects observed at a certain atmospheric seeing are 
blends of several stars. At the same time, typically 
only one star is lensed. In this general case, a de- 
termination of the events' parameters and the analy- 
sis of the detection efficiency of microlensing events 
is very involved and vulnerable to a number of possi- 
ble systematic errors. If the sources are bright, one can 
avoid these problems. First, parameters of the actual 
microlensing events are typically better constrained. 
Second, it is sufficient to estimate detection efficiency 
based on the sampling of the light curve alone. This 
eliminates the need of obtaining deep luminosity func- 
tions across the bulge fields. Red clump giants are 
among the brightest and most numerous stars in the 
bulge, so we focus on these stars here. An optical 
depth determination using all observed microlensing 
events would be desirable, but would require an ac- 
curate calculation of the blending efficiency for non- 
clump stars, demanding much additional input includ- 
ing HST quality images and luminosity functions over 
much of the bulge. 

The structure of this paper is the following. In § 2 
we briefly describe the MACHO experiment as the 
source of the data used here. Section 3 is devoted to 
the selection of microlensing events. In particular, we 
discuss how we select our sample of 62 clump giant 
events (60 unique) and its relation to the catalog of 
over 500 events constructed by Thomas et al. (2004, 
companion paper). In § 4 we test our sample for sig- 



good control over the location of sources they are sensitive to: the 
bulge clump giants dominate over any other possible locations. Dif- 
ference image determinations are sensitive to all sources along the 
line of sight, and thus the determined optical depth is characteristic 
for the direction only and cannot be claimed to be the optical depth 
toward the bulge. A typically employed remedy is to assume that the 
only two populations of sources that matter in star counts are: the 
bulge one and the foreground disk one, and that only bulge sources 
contribute to the optical depth. Following on this assumption, one 
can correct a measured optical depth by a fudge factor that is sup- 
posed to account for the number of inefficient disk sources. It is not 
clear whether such correction is needed and the size of it is very un- 
certain. The corrections applied by both Alcock et al. (2000a) and 
Sumi et al. (2003) increase the optical depth by about 25%. The 
results we quote here do not include these controversial corrections. 
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natures of blending. Optical depth toward the Galactic 
bulge is derived in § 5. In § 6 we consider spatial distri- 
bution of events and discuss the significance of some 
apparent clusters seen in the data. We summarize our 
results in § 7. 

2. Data 

The MACHO Project had full-time use of the 1.27 
meter telescope at Mount Stromlo Observatory, Aus- 
tralia from July 1992 until December 1999. Details of 
the telescope system are given by Hart et al. (1996), 
and details of the camera system by Stubbs et al. 
(1993) and Marshall et al. (1994). Briefly, corrective 
optics and a dichroic were used to give simultaneous 
imaging of a 43' x 43' fields in two bands, using eight 
2048 x 2048 pixel CCD's. A total of 32700 exposures 
were taken in 94 fields (Figure 1) towards the Milky 
Way bulge resulting in around 3 Tbytes of raw im- 
age data and light curves on 50.2 million stars. The 
coverage of fields varies greatly from 12 observations 
of field 106 to 1815 observations of field 119. Note 
that the observing strategy changed several times dur- 
ing the project, so even in a given field the frequency 
of observations changes from year to year. In addi- 
tion, the bulge was not observed at all during the prime 
LMC observing times, so all bulge lightcurves have 
gaps during November-February. 

In this paper, we analyze 7 seasons of Galactic 
bulge data. We select a subsample of 62 clump-giant 
events (60 unique events and 2 duplicate events) from 
337 selection criteria c (to be defined in §3.1) events 
listed in our catalog of more than 500 bulge microlens- 
ing events constructed by Thomas et al. (2004, com- 
panion paper). For the optical depth determination, 
we exclude eleven fields in 300 series (301, 302, 303, 
304, 305, 306, 307, 308, 309, 310, 311). The ex- 
cluded fields are close to the Galactic plane but are 
relatively distant from the Galactic center. As a result, 
they are dominated by the disk stars, which compli- 
cates the morphology of the color-magnitude diagram 
and renders the selection of clump giants much more 
difficult. 

For photometric calibrations, we used global rela- 
tions that express Johnson's V and Kron-Cousins' R 
in terms of the MACHO intrinsic magnitudes 6m and 

r M as: 

V = b M - 0.18(6 M - r M ) + 23.70, (1) 
R = r M + 0.18(6 M - r M ) + 23.41. (2) 
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Fig. 1.— Location of all 94 MACHO fields [includ- 
ing 3 high-longitude fields at (I, b) ~ (18, -2)]. The 
clump regions in fields 301 through 311 are expected 
to be contaminated by disk stars, and thus these fields 
are not included in the optical depth determination. 

For more details see Alcock et al. (1999). The limit- 
ing magnitude in V-band is about 21.5, and a typical 
seeing is 2.2" and almost identical for the blue and red 
filters. 

3. Selection of Events for Optical Depth Determi- 
nation 

3.1. Event Selection 

The microlensing events are selected based on sev- 
eral levels of cuts using statistics calculated for the 
lightcurves. We denote as level-0, all 50.2 million 
lightcurves in our database. Using a small set of vari- 
ability and "bump-finding" statistics, a few percent of 
the lightcurves are advanced to level- 1. The com- 
plete set of over 150 statistics, including non-linear fits 
to microlensing lightcurve shape, is calculated for all 
these lightcurves and around 90000 are advanced to 
level 1.5. Final selection (level-2) is made after fine 
tuning the selection criteria to maximize inclusion of 
genuine microlensing events and minimize inclusion 
of variable stars, supernovae, noise, etc. A more thor- 
ough description of the most useful statistics is given 
in Alcock et al. (2000b). Table 1 gives a brief summary 



3 




10 8 6 4 2 



1 

Fig. 2. — Spatial distribution of events with clump 
giants as sources. Radii of circles are scaled as 
t/e(t, A max ) (See §5). The clump region in fields 301 
through 311 (Fig. 1) are expected to be contaminated 
by disk stars, and thus clump giant events are not se- 
lected in those fields. 

of the selection criteria used here, many of which mea- 
sure signal/noise and goodness-of-fit to a microlensing 
lightcurve shape. 

To differentiate the current set of level-2 cuts from 
the cuts used previously and from the cuts used to se- 
lect LMC microlensing events, we denote the current 
set of cuts as "selection criteria c". It is crucial to use 
exactly this set of selection criteria on the artificial mi- 
crolensing events we create to calculate our detection 
efficiency. 

A total of 337 stars passed selection criteria c, 62 
of which are clump giants (a total of 531 events are 
reported in the companion paper; these were found 
with a less restrictive selection). Figure 2 shows the 
position of the clump events on the sky. Clump gi- 
ant stars are found by final additional cuts on the 
color-magnitude diagram (CMD). This selection dif- 
fers from the selection used in Alcock et al. (1997a). 
A detailed description of the clump giant selection is 
given in the next subsection. The full lightcurves of 
the microlensing events are shown in the companion 
paper (Thomas et al. 2004). 



3.2. Clump Giant Selection 

We pay very close attention to the selection of an 
appropriate clump region. Selection of bulge clump 
giant sources is important since any non-giant or non- 
bulge source stars mistakenly included will distort the 
optical depth result due to being heavily blended or 
being at the wrong distance. Therefore, our goal is to 
select a clean sample of unblended bulge sources. It 
is very hard to distinguish bulge red giants from disk 
clump giants based on the CMD location alone. There- 
fore, we decided to not include bright red stars in our 
clump sample. As a result, we limit clump giants to the 
relatively narrow strip associated with the extinction- 
based selection described below. Also, we introduce 
specific cut to eliminate bright main sequence stars 
that are in the foreground Galactic disk. While clump 
giants form a well defined population the patchy ex- 
tinction towards the Galactic center makes selection 
more complicated. We have attempted to correct for 
this effect in a more sophisticated way than in the past, 
but the possibility for contamination still exists. 

The question, which of our sources may be clump 
giants, is first investigated through the analysis of the 
global properties of the color-magnitude diagram in 
the Galactic bulge. Using the accurately measured ex- 
tinction towards Baade's Window (Stanek 1996 with 
zero point correction according to Gould, Popowski, 
& Terndrup 1998 and Alcock et al. 1998) allows us 
to locate bulge clump giants on the dereddened color- 
magnitude diagram. Such diagram can be then used 
to predict the positions of clump giants on the color- 
apparent magnitude diagram for fields with different 
extinction. Based on Baade's Window data we con- 
clude that unreddened clump giants are present in the 
color range (V - R)q £ (0.40, 0.60), and they con- 
centrate along a line V = 14.35 + 2.0 (V - R) , 
where the zero point of this line is likely uncertain at 
the level of at least 0. 1 mags 19 . We assume that the ac- 
tual clump giants scatter in Vq magnitude around these 
central values, but by not more than 0.6-0.7 mag to- 
ward both fainter and brighter Vq. This range of mag- 

'The (V — R)o color range given here is somewhat bluer than de- 
rived by Popowski et al. (2001a,b) and reported by Popowski et al. 
(2003a). The reason for this change is the detection of anomalous 
extinction law toward the bulge fields probed by microlensing sur- 
veys, which in turn affects the dereddened colors. A possibility of 
smaller than normal Ry coefficient in the bulge direction was sug- 
gested by Popowski (2000) and spectacularly confirmed by Udalski 
(2003) and Sumi (2004) based on the OGLE data. Very similar con- 
clusions can also be drawn from the extinction calibrations derived 
by Popowski, Cook, & Becker (2003b) based on the MACHO data. 
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nitudes is designed to include both the population scat- 
ter and the distance dispersion of clump giants along 
the line of sight. The parallelogram-shaped box in the 
upper left corner of Figure 3 is approximately defined 
by the above assumptions. To be exact, we define this 
parallelogram to have the following vertex points in 
[(V - R)o,V ] space: (0.40, 14.52), (0.40, 15.88), 
(0.6, 16.28), and (0.6, 14.92). With the assumption 
that the clump populations in the whole bulge have 
the same properties as the one in the Baade's Window, 
the parallelogram described above can be shifted by 
the reddening vector to mark the expected locations of 
clump giants in different fields. The solid green lines 
are the boundaries of the region where one could find 
the clump giants in fields with different extinctions. 

A universal selection of clump giants outlined 
above (see also Popowski et al. 2001a, 2003a) allows 
some contamination by main sequence stars in high 
extinction areas of the sky. Popowski et al. (2001a,b) 
tried to remedy this problem introducing a universal 
color cut (V — R) > 0.7. However, the stars in 
high-extinction regions may become so red that the 
main sequence enters so-defined clump region from 
the blue. On the other hand, moving this generally 
fixed limit of (V — R) = 0.7 to the red would re- 
move legitimate clump giants from low-extinction re- 
gions. Therefore, one is driven toward introducing a 
non-universal clump selection that makes use of the 
properties of the CMD at a given sky location. 

In our non-universal selection, we decided to keep 
our previous general procedure of selecting clump gi- 
ants intact except for 1) the adjustments coming from 
the different form of the reddening-law adopted here, 
2) addition of a blue color cut tied to the characteristics 
of the CMD. We checked that the region that is small 
enough to be sensitive to the local changes of extinc- 
tion and large enough to form a reliable CMD approx- 
imately coincides with the, so-called, tile in the MA- 
CHO database. Each tile is a 4 by 4 arc-minutes patch 
and contains a few thousand identified stars. As our 
bulge fields contain more than 10000 tiles, it would be 
very impractical to define blue edges of clump regions 
on an individual basis. Therefore, we formed a training 
set of 240 tiles (in places corresponding to suspected 
microlensing candidates from the first 5 years of the 
experiment, but rejecting all tiles in 300 series fields 
and keeping only one copy of repeated tiles). This 
set of tiles is fairly representative and contains regions 
with a range of stellar densities and stellar extinctions. 
We decided to set a blue clump limit at (V — R) half 



the way between the central clump and main sequence 
over-density at the V-magnitude of the clump. The 
240 CMDs were visually inspected and the most likely 
clump boundary was selected. The next task was to 
relate these boundaries to the global properties of indi- 
vidual CMDs. We computed first four moments of the 
(V— R) color of each CMD (no luminosity weighting). 
We used the mean and dispersion as the only parame- 
ters in the linear fit. Mathematically, we requested: 

(V - i?)boundary = OL + (3 (V - R) CMD 

+ 1 f{a{V-R)cMD), (3) 

where we tested three types of f(x), namely f(x) = x, 
f(x) = x 2 , and f(x) = ln(x). All forms of equa- 
tion (3) produced very similar results, and we adopted 
f(x) = ln(x), requesting that the clump giants be red- 
der than: 

(V - i?)boundary = (-0.001 ± 0.033) 

+ (0.872 ±0.014) (V-R) CMD 
+ (-0.042 ± 0.012) In {a(V - R)cmd) ■ (4) 

The errors reported in equation (4) were normalized 
to produce x 2 /d.o.f. = I. 20 The clump blue bound- 
ary determined by visual inspection that deviates most 
from the fit reported in equation (4) is at 3.28ct away. 
In general, the number of 3er points is not very differ- 
ent than expected from a normal distribution of errors, 
and therefore we reject no points as outliers. More- 
over, in the cases with largest deviations, the blue lim- 
its suggested by equation (4) are equally reasonable as 
originally selected boundaries, and, therefore, do not 
hint at any small subset of CMDs that do not obey a 
general relation. The average scatter of the visually- 
selected (V — i?)boundary around relation (4) is below 
0.02 mag, which, being an order of magnitude smaller 
than a typical range of clump colors, is more than sat- 
isfactory. 

Our final set of clump selection criteria are de- 
scribed by the following equations: 

Hasc > 15.0 and Vb asc < 20.5, (5) 

V hase > 4.2 (V - i?) basc + 12.4, (6) 
V base < 4.2 (V - i?) basc + 14.2, (7) 

'The correlation matrix is: 

1.0000 -0.9084 0.9720 \ 

-0.9084 1.0000 -0.7872 

0.9720 -0.7872 1.0000 / 
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boundary! (8) 

field < 301 or field > 311, (9) 

where subscript "base" referring to V^-magnitude and 
(V — i?)-color indicates baseline magnitudes, i.e. the 
ones in the limit of no microlensing-induced amplifi- 
cation. The cut on bright sources in (5) is intended 
to eliminate events that are in the foreground or may 
have saturated photometry, the cut on faint sources is 
designed to avoid stars with uncertain photometry or 
highly affected by blending. In practice, the cut on 
bright sources from (5) is always weaker than the com- 
bination of cuts (6) and (8). Cuts (6) and (7) deter- 
mine the reddening channel that marks the possible 
locations of clump giants. Cut (8) admits only stars 
with red enough color, which removes contamination 
by foreground main sequence stars. Finally, cut (9) ex- 
cludes 1 1 disk dominated/contaminated fields, namely 
the ones in the 300 series. 

The slope of relations (6) and (7) has the interpreta- 
tion of the coefficient of selective extinction Rv,vr = 
Ay /E{V — R). The coefficient we use here Rv,vr = 
4.2 was found by Popowski et al. (2003b) in Baade's 
Window. We checked that clump selection based on 
Rv,vr = 4.5 produces very similar results. On the 
other hand, Rv,vr — 5.0 seems to be too steep. De- 
spite this, Popowski et al. (2001a) selection was quite 
successful, because their intrinsic clump color was as- 
sumed to be redder and they had very few events in 
high-extinction regions, which are primarily affected 
by this change. Finally, let us note that the properties 
of extinction law and the zero-point needed to deter- 
mine absolute extinction in the bulge are not known 
precisely, which could have negatively affected our se- 
lection. This was not the case. We visually checked in 
the number of regions with different level of redden- 
ing that the operational procedure defined in equations 
(5)-(9) properly picks up clump giants in a wide range 
of color-magnitude diagrams. 

Clump region defined by conditions (5)-(9) is plot- 
ted in Figure 3 and surrounded by bold green lines. 
In addition to the boundary of clump region, we also 
indicate cut (5) by blue lines, cuts (6) and (7) by red 
lines, and cut (8) by dotted magenta line, with arrows 
indicating its non-universal, CMD-specific character. 
The clump events are marked as filled triangles. The 
events inside bold green lines marked as open triangles 
belong to fields in 300 series, and are excluded by cut 
(9). Alcock et al. (1997a) clump region is presented 
for comparison. 
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Fig. 3. — Color-magnitude diagram of MACHO ob- 
jects. Triangles are events and filled triangles are 
clump events. The underlying population is shown in 
yellow. 

The properties of the 62 lensed clump giant stars 
as well as the microlensing fit parameters are given in 
Tables 2 and 3. In these and the following tables V 
and (V — R) will refer to baseline magnitudes - we 
drop subscript "base" to avoid clutter. We note that 
two events are marked with letters ('a' or 'b') indicat- 
ing they are duplicates of another event (duplicate also 
marked 'a' or 'b'). In general, this is caused either by 
the same physical event being detected in an overlap- 
ping field, or by a nearby star receiving flux from the 
actual event due to incorrect flux sharing in the pho- 
tometric PSF fitting code. Both duplicate events listed 
in Tables 2 and 3 result from field overlaps. For the 
optical depth calculation we will count both of the du- 
plicate events since field overlaps increase both the to- 
tal number of stars monitored and the number of events 
detected proportionally, and since the field overlap and 
number of duplicates is fairly small. Finally we note 
that 5 of the clump giant microlensing events are po- 
tential binary events, marked with a * symbol. 

In Figure 4 we graphically summarize the basic 
properties of clump events. Most of the panels are 
self-explanatory. Let us only remark about the distri- 
bution of the number of events versus the time of max- 
imum amplification to- The peaks in the histograms 
show 7 observing seasons analyzed here. The number 
of detected events varies mostly due to Poisson noise 
with the noticeable exceptions of the second and last 
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Fig. 4. — Basic properties of the clump sample. 

seasons. The numbers of events are particularly low 
there due to an observing strategy that targeted very 
few fields in search for short-duration events. 

In Figure 5 we test for possible systematic effect 
in microlensing parameters. In particular, in the upper 
four panels we plot impact parameters and Einstein di- 
ameter crossing times as a function of Galactic coor- 
dinates. One sees no dependence of parameters upon 
position. The left lower panel suggests that there is 
no correlation between the impact parameter and du- 
ration, and the right lower panel shows that the impact 
parameter is not related to the event color. 

Finally, we note that many of the events selected 
here were alerted on by us earlier in the experiment. In 
addition, the MACHO fields overlap with fields from 
other microlensing experiments. We cross-referenced 
the clump sample with our own alert events and with 
events from EROS, MOA, and OGLE experiments. 
We found no counterparts in the MOA set (as expected 
from non-overlapping observing seasons). Alterna- 
tive designations from the MACHO Alert system and 



Fig. 5. — Microlensing properties of events as a func- 
tion of position, duration, or color. 

counterparts from EROS and OGLE are listed in Ta- 
ble 4. 

4. Is Blending Negligible for the Clump Sample? 

We make an analytic estimate of the effect of blend- 
ing on our clump analysis using the scaled luminos- 
ity function for Baade's Window from Holtzman et al. 
(1998). We find that even under very conservative as- 
sumptions, the number of highly blended clump events 
should not exceed 25%. Moreover, this level of con- 
tamination affects the optical depth at the level not ex- 
ceeding 10%, and for some photometric code behav- 
iors the contamination may be completely negligible. 
The details are discussed in Appendix A. To check 
whether the sample of selected clump giants meets our 
expectations of being weakly affected by blending we 
applied several tests. 

First, we reviewed the color light curves of all 
events with particular emphasis on the peak region. 
For clump giants, blends with bluer stars are more 
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likely than the blends with stars of identical color 
and so this test should be sensitive to the major- 
ity of potential contaminants The light curves of all 
clump events showing blue and red magnitudes and 
the difference between them for simultaneous ob- 
servations are shown in Figure 6. Here we display 
two example events, whereas the entire sample is 
available in the electronic version and will be avail- 
able on the World Wide Web upon the acceptance of 
this paper 21 . The red lines superposed on the data 
are the best fits obtained under the assumptions of 
no blending, and in most cases they properly rep- 
resent observational points. We examine the light 
curves for any deviations from the unblended fits. 
We note that events 109.20640.360, 113.19192.365, 
176.18826.909, 401.47994.1182, 401.48052.861, 
401.48167.1934, 402.47856.561, and 403.47491.770 
either lack or have extremely sparse color coverage in 
the peak region. Therefore, any blending-related in- 
formation that can be extracted from their light curves 
is very limited. The events for which there is at least 
a slight possibility of achromatic signal are listed in 
Table 5. Several events have deviation from fitted light 
curves that are due to their exotic character: either bi- 
narity or parallax effect. Only 1 1 out of 62 events are 
not immediately explained by known effects. In addi- 
tion, 3 of them have either asymmetric or caustic type 
signal which is more likely an indication of their bi- 
nary character than blending. Therefore, only 8 events 
(13%) can be suspected of significant blending based 
on the visual inspection of their light curves. 

Next we examined the color close to the center of 
the peak (t — to < i/ y/2) compared to the baseline of 
the lightcurve in a more quantitative fashion. We start 
with a region of width t/8 centered on to and expand 
it until we have at least 2 points with simultaneous ob- 
servations in each filter. Using these points in each 
lightcurve we calculate the difference in flux between 
the peak and the baseline in each filter (AV or AR). In 
Figure 7 we plot the ratio AV/ AR versus A max . For 
blending-free microlensing events AV/ AR should be 
unity since gravitational lensing is wavelength inde- 
pendent. Of the 53 events that had simultaneous ob- 
servations in both filters in the peak, 20.7% were more 
than lcr away from unity. This is consistent within er- 
rors with no blending in our sample. There are two 
events (108.18951.593 and 403.47610.576) which dif- 
fer by more than 6 sigma from unity, which might be 
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indicative of blending. One (108.18951.593) is likely 
blended with a faint star of a very different color, 
and examination of the lightcurve of the second one 
shows that it has many points far off the microlensing 
fit (403.47610.576), probably indicating intrinsic vari- 
ability or problems with the photometry of this star. If 
we ignore these two events we have 17.6% more than 
lcr from unity and the distribution is roughly Gaussian. 

As our third test, in Figure 8, we plot the cumula- 
tive distribution of the impact parameters u m i n and a 
uniform distribution (straight line) of u m j n between 
and 0.826. This distribution is what is predicted for 
microlensing events with magnifications between in- 
finity and 1 .5 as imposed by our selection criteria. We 
find good agreement, with the Kolmogorov-Smirnov 
D = 0.107 for 55 events (60 unique events minus 5 
binaries). Under the hypothesis that all events are mi- 
crolensing, this gives a probability of 53% of obtaining 
a value of D this large or larger. For this analysis we 
made a correction for the efficiency of each event, but 
even without making this correction we get excellent 
agreement (D = 0.103 for a probability of 58%). This 
suggests that u m i n values from standard microlensing 
fits are not seriously biased by blending. 

Finally, we perform so-called blend fits. That is we 
fit every flux curve with a formula 

F(t) = F hasc [(1 - /) + fA(t)] , (10) 

where A(t) is the amplification of standard point lens 
microlensing, fbase is the baseline flux of an object 
unaffected by the microlensing magnification, and / 
is the blending fraction. For a completely unblended 
source one expects / = 1. We fit all curves using 
MINUIT (CERN Lib. 2003), and report the results 
together with the parabolic errors in Table 6. In our 
fitting setup, the blending fractions for an event can 
take different values in blue and red filters. The up- 
per panel of Figure 9 shows the distribution of blend 
fractions for the events from the CGR (a group of 9 
fields used for the determination of the average opti- 
cal depth, which is defined in §5). As expected, many 
events have blend fractions close to 1.0, i.e. consis- 
tent with no blending. At the same time, however, the 
results look quite alarming. The distribution of blend 
fractions is very broad and there is a number of events 
with small blend fractions that indicate heavy blend- 
ing. To understand the meaning of this result we per- 
form Monte Carlo simulations of the expected distri- 
bution of recovered blend fractions for completely un- 
blended clump events from a sample identical to the 
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one presented in the upper panel of Figure 9. We cre- 
ate light curves of unblended events with durations and 
Mmin values identical to the ones from the real sample. 
We produce a few hundred light curve realizations for 
each event and compute for them blend fit parameters. 
Thus for each real event we obtain an expected distri- 
bution of recovered blend fractions if the original event 
was unblended. We sum those distributions together 
to obtain the expected blend fraction distribution for 
the real sample. We plot the resultant distribution in 
the lower panel of Figure 9 (solid line). Therefore, 
even if there is no blending in our clump sample, the 
blend fit procedure will predict many blended events. 
This problem has been studied before and is a result 
of sampling and data quality (Wozniak and Paczyhski 
1997). The shapes of the distributions from Figure 9 
are similar, which suggests that our null hypothesis 
that the clump sample is unblended may be correct. 
For comparison, we also plot the recovered blend frac- 
tions for the sample with input blend fractions of 0.5 
(dashed line). It is clear that the peak follows the true 
blend fraction, which again suggests that our sample is 
mostly unblended. In addition, the distribution result- 
ing from the input blend fractions of 0.5 is narrower 
than the true one and lacks extended wing toward large 
/ values. However, we also note that the real sample 
seems to have some excess of events with very small 
blend fractions, inconsistent with neither of the two 
distributions shown in the lower panel. 

To further test possible contamination of the clump 
sample by highly blended events, we notice that the 
blending fractions themselves do not tell the entire 
story. There is no reason to reject a non-blending hy- 
pothesis as long as the blend fractions are consistent 
with f = 1 . Therefore, we compute the deviation 8 of 
a blend fraction from f = 1 in units of the error in the 
blend fraction: 

«-<^±2, (ii) 

where 8 assumes negative values for / < 1.0 and posi- 
tive values for / > 1.0. 22 It is obvious from the defini- 
tion that 8 will not significantly differ from for events 
consistent with no blending. The events with large de- 



The blend fractions with / > 1.0 seem unphysical, but even in 
the case of negligible observational errors they can result from local 
fluctuations in the sky level. Such fluctuations will typically pro- 
duce / < 1.2. However, in the presence of noise / > 1.2 solutions 

should not be rejected as long as they are consistent with physically 
motivated / values. All solutions in our sample meet this basic con- 
dition. 



viations can be classified as suspicious. We plot the 
distributions of <5bi U c for the real and simulated sam- 
ples in Figure 10. Again, the distributions are not very 
different from each other, but we cannot exclude the 
possibility of significant blending for some events in 
our original sample of 60 unique clump events. We 
conclude that, given the quality of our data, the blend 
fit parameters we find for events with a single lens can- 
not be used to make a definitive determination of the 
amount of blending in our sample. Nonetheless, since 
our ultimate goal is to obtain a reliable estimate of the 
optical depth, we use blend fit results to select a subset 
of our sample, which we design to be fully consistent 
with no blending. This extremely conservative subset 
consists only of the events with |<5biuc| < 2.0. It is en- 
couraging that this subsample contains very few events 
from Table 5. The analysis of this verification sample 
conducted in the next section will provide an essential 
check on the systematic error in our primary optical 
depth determination based on the entire clump sample 
assumed to be blending-free. 

Finally, we notice that blending fits to caus- 
tic crossing events with binary lenses are substan- 
tially better constrained than the fits to events with 
single lenses. Our clump sample contains 4 bi- 
naries that have been previously analyzed by Al- 
cock et al. (2000c): 108.18951.593 (97-BLG-28), 
118.18141.731 (94-BLG-4), 40 1 .48408 .649 (98-BLG- 
14), and 402.47864.1576 (97-BLG-41). The first three 
have fits in Alcock et al. (2000c, Table 3), and 97- 
BLG-41 has two possible solutions: binary with a 
planet suggested by Bennett et al. (1999) and rotat- 
ing binary discussed by Albrow et al. (2000). Out of 
16 bulge binaries that have blend fits in Alcock et al. 
(2000c) only 6 are fully consistent with no blending 
(94-BLG-4, 96-BLG-4, 97-BLG-l, 97-BLG-28, 98- 
BLG-14 ("dashed solution"), and 108-E). On the other 
hand, the four clump binaries for which we have fits 
have blend fractions very close to 1 and are consistent 
with no blending (see their blend fractions in Table 7). 
The probability that the clump binaries share the same 
blending characteristics as the entire sample of bina- 
ries is P(4) ee (^)(6/16) 4 (10/16)° = 0.020. There- 
fore, we have a 2.3 a indication that our clump sample 
may be affected by blending to negligible extent. 

5. The Optical Depth 

We propose to use an estimator of the optical depth 
which is different from, but closely related to, the one 
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we used previously. The estimator used in our previous 
papers is: 

■^events 7 

where £ is time to cross the Einstein ring diameter, E 
is the total exposure in star-days (equal to the number 
of stars times time span of observations), and e(t) is an 
efficiency for detecting an event with a given t. 

A formula which gives the same answer in the limit 
of large numbers of events, but which we find more 
computationally efficient is: 

-i N events > 

t= e £ 7(t- a r (13) 

- n C V z> rnax,z/ 
i— 1 

where ti is the measured time for which the magnifi- 
cation is above A = 1.34, that is the time the line- 
of-sight is within the Einstein ring radius. Appendix 
B gives the derivation of both of these formulae and 
makes the case for the use of equation (13). Note that 
due to our cut A max > 1.5, we must multiply the value 
derived from equation (13) by a correction factor of 
1.09. This is explained in Appendix B. 

There are five likely binary events in our sample. 
Not only is the parameter determination more difficult 
in their case, but our detection efficiency pipeline de- 
scribed below is designed for the events with a single 
lens. Therefore, it is not immediately clear how to in- 
clude binaries in the analysis. They cannot be simply 
omitted because this would lead to an underestimate of 
the optical depth. For simplicity we used for binaries U 
as found from single lens fit. We found this quite accu- 
rately measured the time above A = 1.34, and agrees 
quite well with what is found from the binary lens fits 
(Alcock et al. 2000c) 23 . Even if make a combined sys- 
tematic error of 50% in the estimate of the time above 
A = 1.34 and detection efficiency for binary events, 
this will result in only about 4% error in the optical 
depth because binary events constitute only 1/12 of the 
entire sample. Another problem could occur if our cuts 
selected against binaries due to non-standard shapes of 
their light curves. Taking the fraction of binaries for all 
events (selected with much more relaxed cut) from the 
companion paper, we find that the expectation number 
of binaries in the clump sample is six. We see five. If 
we count one missing binary as a source of bias in the 



We notice that the values of t found from the single lens fits are not 
very useful for these events. 



optical depth, then we predict the optical depth under- 
estimate to be at the level of (1/60) x 100% = 1.5%. 
Both effects are small compared to the statistical error 
in the optical depth. 

The sampling efficiencies were obtained with the 
pipeline that has been previously applied to the LMC 
data (for a description see Alcock et al. 2000b). In 
brief, a random subset of 1% of all lightcurves was 
selected and artificial microlensing light curves with 
different parameters were added to these lightcurves. 
Then the same analysis used to select real events was 
applied to this set. Efficiencies were computed as the 
number of events recovered divided by the number 
of events input as a function of both the input dura- 
tion and the input maximum magnification of the in- 
put event. Since the number of clump giants in each 
1% sample field is not large, several passes were made 
through each 1 % data base, until enough statistics were 
gathered. We averaged over A max to report efficien- 
cies for each field as a function the event duration (i), 
and we also did special runs using the measured t j and 
^4max,i for each real clump giant microlensing event. 
This last analysis gave us e(U, j4 maXl i) for use with ti 
in equation (13) above. Figure 1 1 shows the results of 
our efficiency calculations for the fields with events, 
and Table 3 shows the efficiency calculated for each 
event. Relative sampling efficiencies for all 94 bulge 
fields are given in Table 8. The results in fields in the 
300 series are not used in the optical depth analysis 
and they may be affected to higher extent by system- 
atic effects. 

Using the above method and the clump events from 
Table 3, we present in Table 9 optical depths for 
each field that contains clump giant events. The cen- 
ter position of each field, the number of clump giant 
source stars (in the text designated shortly by Nf), and 
the number of recovered microlensing events are also 
given. The number of days, T b s = 2530, used to find 
each field's exposure E = NfT ^ s was set by the ef- 
ficiency calculation 24 . Since the number of events in 
each field is so small the error in the optical depth is 
dominated by Poisson statistics. We compute these er- 
rors using the formula given by Han & Gould (1995). 

We show in Figure 12 the spatial distribution of op- 
tical depths on the sky. In the left panel we give the 
optical depth value at the center of each field with 



The number of days T b s is equal to the sum of the allowed to span 
(2850cZ - 419d = 2431d) and two ~ 50-day buffer regions added 
at the beginning and the end of the observing interval. 



10 



clump events, and in the right panel we present the 
color representation of the same results and indicate 
the number of detected events. We note the anoma- 
lous character of field 104, with an optical depth of 
(8.8 ± 3.7) x 10~ 6 , almost two sigma above the av- 
erage. The results for individual fields are very un- 
certain but taken as an ensemble they suggest that the 
optical depth has a substantial gradient in Galactic lat- 
itude. This property is more apparent in Figure 13, 
where we binned the data on a 0.5 deg scale in the 
fe-direction and 1.0 deg scale in the /-direction. We 
then fit straight lines to these binned data. Based on 57 
clump giant events selected from observations of the 
1260000 clump giant source stars within 5.5 degrees of 
the Galactic center we find the optical depth gradients 
of (1.06 ± 0.71) x 10 _6 dcg _1 in the 6-direction and 
(0.29±0.43) x lO^dcg -1 in the /-direction. It is clear 
from the errors that the gradient in the /-direction is 
consistent with and with different small values close 
to (both positive and negative). 

Many fields do not contain any clump events. In 
those cases we derive an upper limit on the optical 
depth from 

lim.CL = CL [ n Y^=l mt 'iihSi) \ 

(14) 

where e/ is the efficiency in field /, ej is the effi- 
ciency for i-th event in its field of origin, and ^ CL is 
the multiplication factor dependent on the confidence 
level. Equation (14) is derived in Appendix C. We use 
v lcr = 1.8379 when we compute the upper limits on 
the optical depth given in Table 10. Those upper lim- 
its were computed without any direct reference to the 
Galactic model. However, they rely on the assump- 
tion that the intrinsic duration distribution in empty- 
fields is identical to the one in fields with events, and 
this approximation can break down for some Galactic 
models. In addition, let us note that the limits from 
Table 10 are rather conservative in the sense that the 
true values may be lower. This is caused by the fact 
that we used durations of all events when we applied 
equation (14). We show in the next section that field 
104 tends to have longer events. If events in field 104 
were excluded from the summation in equation (14), 
then the limits would go down a little. 

Since the number of events in individual fields is 
typically very small, microlensing studies customar- 
ily give average optical depth based on many fields 
to beat the Poisson noise. The main problem with 



this approach is the fact that the interpretation of the 
average optical depth and position is difficult due to 
non-linear variations of r over the area of the survey. 
Decreasing the size of the region for averaging results 
in reappearance of the Poisson noise issue. Here we 
deal with these competing effects by selecting a very 
compact region that contains 2/3 of all clump events. 
We call this region Central Galactic Region (CGR) to 
stress its close proximity to the Galactic center (but 
we would like to clarify that there are no MACHO 
fields that include the Galactic center itself). CGR 
covers about 4.5 sq. degrees and is composed of 9 
fields: 108, 109, 113, 114, 118, 119, 401, 402, 403. 
Using 42 events in CGR, we find an average opti- 
cal depth of t cgr = 2.17±^;|| x 10~ 6 at (l,b) = 
(1?50, — 2?68) 25 . For this set of fields, an exposure is 
E = 1.869 x 10 9 star-days (2530 days times 738850 
clump giant stars). In this case the error was found by 
Monte Carlo simulations as described by Alcock et al. 
(1997b), and is somewhat larger than the simple Pois- 
son error (which would be ±0.34 x 10~ 6 ), but consis- 
tent with the formula of Han and Gould (1995), which 
gives ±0.42 x 10~ 6 . We note that the result found 
using equation (12) and an interpolated efficiency as 
a function of t is nearly identical (2.16+o' 3 g x 10~ 6 ). 
Field 104, with the largest optical depth, falls just out- 
side of the CGR. To test how sensitive the results are to 
the detailed shape of the region selected for the deter- 
mination of the average optical depth, we create a rect- 
angular region "CGR+3", which is composed of CGR 
plus the three adjacent fields at the same galactic lati- 
tudes (fields: 176, 104, 105). The optical depth in such 
an extended region is t C gr+3 = 2. 37toi39 x 10~ 6 at 
(l,b) = (1?84, -2? 73), entirely consistent with our 
principal result. Also, we note for completeness that 
when we use all 2.44 million clump giants and all 62 
events the optical depth over all 83 analyzed fields is: 
r a n = (1.21 ±0.21) x 10~ 6 at (l,b) = (3?18, -4?30), 
but this number is hard to interpret given non-linear 
(and often rapid) variations of the optical depth with 
Galactic position mentioned earlier. 

Finally, we return to the question of blending and its 
effect on the optical depth. To test the sensitivity of the 
optical depth to possible systematic errors, we analyze 
the verification sample of events. To avoid possible 
systematic errors we first exclude all binaries. Then 
we keep only the events that have |<5bi uo | < 2.0, which 

The unweighted average position for 9 CGR fields is (I, b) = 

(1?55, -2? 82) 
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makes them very consistent with no blending. We list 
the events that pass this additional cut in Table 1 1 . We 
list events that belong to CGR on the left and the events 
in the other fields on the right. 

The analysis on this sample proceeds in a fashion 
that is similar to the main analysis but has a noticeable 
difference. Our verification sample is constructed by 
utilizing the cut that uses blend fit results. Applying 
the same procedure to our massive efficiency simula- 
tions is beyond the scope of this paper. Therefore, we 
use the previously determined efficiencies and the fol- 
lowing estimator of the optical depth: 



r vcr events 

- N vcr 

events 
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i event; 
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(15) 



where N°"^ ts and N^ nts are the numbers of events 
in a given field or region counted in the original and 
verification sample, respectively. In designing the op- 
tical depth estimator given in equation (15) we as- 
sumed that blending (or other systematic effect we 
would like to test) does not change the number of de- 
tected events, but it may substantially affect the event 
parameters and may render the sampling efficiencies 
inappropriate. However, both the parameters and ef- 
ficiencies can still be used for the verification sample. 
The factor N°"^ ts / N^ nts is introduced to compen- 
sate for the fact that we remove some events without 
replacing their optical depth contributions with any- 
thing else. 

The optical depth values in individual fields de- 
rived using equation (15) for our verification sample 
are given in Table 12. The first four columns give 
the field number, Agents' tne optical depth derived 
using equation (15), and its error. For comparison, 
columns 5 and 6 give N°^ tB and the optical depth 
based on equation (13) derived earlier. The results 
from the verification sample agree very well with our 
previous estimates. As the last step we limit the ver- 
ification sample to the CGR 26 . The average optical 
depth based on this verification sample of 22 events 
is t<5g R = 2 -42±o.58 x 1(T 6 , in excellent agreement 
with our main result. 



26 It may seem surprising that only 59% or 22 events out of 37 non- 
binary events in the CGR pass |<5bi ue | < 2.0 cut. We note how- 
ever that the corresponding ratio for our simulated unblended sam- 
ple from the lower panel of Figure 10 is 80%. This difference is at 
least partly explained by the fact that the Monte Carlo simulations 
do not include exotic events (e.g., parallax events) which are present 
in the real sample and can cause the blend fits to falsely indicate 
blending. 



In summary, we see no evidence that our optical 
depth results are strongly affected by blending. There 
are two likely explanations of why this is the case. 
First, it is possible that the clump sample is to a large 
extent unaffected by blending. Second, substantial 
blending is a result of a large number of unresolved 
sources. In such situation, the bias in the recovered 
event duration is countered by the bias in the inferred 
number of sources. At the moment, we do not have 
enough information to distinguish between these two 
possibilities, but our optical depth determination ben- 
efits from their fortunate properties. 

6. Clustering of Events and Concentration of 
Long Events in Field 104 

The optical depth for field 104 is larger than any 
other field and we would like to investigate whether 
this is a statistical fluctuation or an indication of some- 
thing unusual about this region of the Galaxy. The 
large optical depth is related to the fact that 9 of the 62 
clump giant events are in field 104, and because there 
is also a high concentration of long-duration events in 
this field. Four out of the 10 events longer than t = 100 
days are in field 104, including the 2 of the 3 longest. 

We investigate how statistically significant this con- 
centration is. We do not account for the change in 
the detection efficiency of events with different dura- 
tions in different fields, so we place only a lower limit 
on significance of the difference between duration dis- 
tributions in field 104 and the others. The efficiency 
for detecting long events is similar in most fields with 
events observed for 7 years, because this does not de- 
pend strongly on the sampling pattern. The detection 
of short events will be lower in a sparsely sampled 
fields. Therefore, the number of short events in some 
of the fields used for comparison may be relatively too 
small with respect to a frequently-sampled field 104, 
but this is only going to lower the significance of the t 
distribution difference computed here. In conclusion, 
the analysis of event durations uncorrected for effi- 
ciencies should provide a lower limit on the difference 
between field 104 and all the remaining clump giant 
fields. 

For the significance test, we use the Wilcoxon's 
number-of-element-inversions statistic. First, we di- 
vide events into two populations: events in field 104 
and all the remaining ones. We test two types of pos- 
sible separations into those two populations: A. we se- 
lect all events in field 104 to the first sample (9 events) 
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and all the remaining events except duplicates from 
field 109 to the second sample (51 events), B. we se- 
lect all events in fields 104 and not in any other field to 
the first sample (7 events) and all the events in the re- 
maining fields the second sample (53 events). In both 
cases the entire sample of 60 unique clump giants can 
be recovered as the sum of populations in 104 and not 
in 104. 

Case A and B are then each tested separately ac- 
cording to the following procedure. We order the 
events in the combined sample of two populations 
from the shortest to the longest. Then we count how 
many times one would have to exchange the events 
from field 104 with the others to have all the 104 field 
events at the beginning of the list. If N\ and N 2 desig- 
nate numbers of elements in the first and second sam- 
ple, respectively, then for Ni > 4, N 2 > 4, and (Ni + 
N 2 ) > 20, the Wilcoxon's statistic is approximately 
Gaussian distributed with an average of NxN 2 /2 and 
a dispersion a of yjN x N 2 {N[ + N 2 + 1)/12. 

In case A, the Wilcoxon's statistic is equal to 324, 
whereas the expected number is 230 with an error 
of about 48. Therefore the events in 104 differ (are 
longer) by 1.96<7 from the other fields. In case B, the 
Wilcoxon's statistic is equal to 282, whereas the ex- 
pected number is 1 86 with an error of about 43 . There- 
fore the events in 104 differ (are longer) by 2.22er from 
the other fields. Both divisions of the clump sample 
into field 104 and the remaining fields lead to the con- 
clusion that events in fields 104 are inconsistent at the 
2 a level with being drawn from the same parent pop- 
ulation as events in the remaining fields. We obtain 
similar discrepancy between the time scales of events 
in field 104 and the remaining fields when we use all 
318 unique events from selection criteria c listed in Ta- 
ble 3 of the companion paper. This result is somewhat 
less significant than the one derived by Popowski et 
al. (2001a). We note, however, that similar analysis 
of an almost independent Alcock et al. (2000a) DIA 
sample also suggests anomalous character of duration 
distribution in field 104 (Popowski 2002). We con- 
clude that the unusual character of field 104 should be 
investigated with additional observations but we can- 
not completely exclude the possibility that this effect 
is due to a statistical fluctuation. 

In addition, there seem to be clusters of events (of 
all durations) on the sky in fields 104, 108, and 402. 
If microlensing events are clustered on the sky above 
random chance it has important consequences. It could 
indicate clustering of lenses, perhaps in some bound 



Galactic substructure or very young star forming re- 
gion. The long duration of the events could indicate 
that the substructure contained heavier objects such as 
black holes. Alternatively, it could indicate a concen- 
tration of orbits through the Milky Way and a place 
where the orbital speeds were slower than normal. 

In order to test the significance of the apparent 
clustering of microlensing events we devised a Monte 
Carlo test in which we simulated 10000 microlens- 
ing experiments each detecting 60 unique clump gi- 
ant events. In each experiment we found densest clus- 
ter of 3 events and the densest cluster of 4 events and 
recorded the diameter of the cluster. We then com- 
pared the diameters of the densest 3-cluster and dens- 
est 4-cluster in the actual data to the density distribu- 
tions formed by our Monte Carlo experiments. This 
gives us the probability of finding by chance a 3-cluster 
(or 4-cluster) as small as we found in the actual data. 
In both the data and in the simulations we remove any 
star that is found two or more times (e.g., in an over- 
lap region), since these give a false measure of dense 
clustering. 

In performing the simulations we selected 60 clump 
giant stars at random from the 1 percent data base so 
that the spatial distribution of lensing sources matched 
that in the data. We also weighted each event by the 
50-day efficiency of that field (Table 8), so that fields 
with low sampling efficiencies were properly under- 
represented. Finally we considered various gradients 
in optical depth across our fields. As one limit we use a 
uniform optical depth across our fields and as the other 
we use a steep linear gradient which results in the op- 
tical depth changing from a maximum at the Galactic 
center to at \b\ = 4°. This later case, is roughly con- 
sistent with, but somewhat steeper than, the gradient 
shown in Figure 13. The actual gradient should be be- 
tween these two extremes. The second case increases 
the chance of randomly finding a cluster with respect 
to the first case. 

With no optical depth gradient we find a probabil- 
ity of 1.9% of finding a 3-cluster as dense as in the 
data, and 4.6% of finding a 4-cluster this dense. For 
the steep optical depth gradient, we find probabilities 
of 8.7% and 27% respectively. Unfortunately these re- 
sults are not very conclusive. The clustering of mi- 
crolensing events we find in our data is marginally sig- 
nificant, but certainly not compelling. We note for 
completeness that the densest clump 3-cluster in real 
sample had a diameter of 0.0453 degrees and the dens- 
est 4-cluster had a diameter of 0. 1 1 degrees. 
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We also performed the clustering analysis on the 
entire sample of 318 unique microlensing events se- 
lected with criteria "c". Here the densest 3- and 4- 
clusters in real data had the diameters of 0.0272 and 
0.0497 degrees, respectively. Note that for random sur- 
face distribution of points the size of the smallest clus- 
ter scales with the number of events N as iV -1 / 2 . We 
conclude that the above diameters are very consistent 
with the clump only result since they produce rescaled 
diameters of 0.0626 and 0.114 for 3- and 4-clusters re- 
spectively when multiplied by ^/318/60 = 2.302. For 
the entire sample of 318 unique event and for uniform 
optical depth, the probability of finding a 3-cluster as 
dense as in the data is 7.3% and the probability of find- 
ing a 4-cluster this dense is 3.6%. For the steep optical 
depth gradient the probabilities are 36% and 32% re- 
spectively. So for the entire sample of microlensing 
events we also find no strong evidence of clustering on 
the sky. 

7. Conclusions and Discussion 

Using 7 years of MACHO survey data, we pre- 
sented a new determination of the optical depth to mi- 
crolensing towards the Galactic bulge. We selected 
the sample of 62 microlensing events (60 unique) on 
clump giant sources (out of more than 500 total) and 
performed a detailed efficiency analysis. Using a sub- 
sample of 42 clump events detected during monitor- 
ing of about 739000 clump giant sources concentrated 
in just 4.5 deg 2 , we found r = 2.17±g;|l x 10~ 6 at 
(I, b) = (1?50, -2?68), consistent with theoretical ex- 
pectations. 

We can make a comparison with previous MA- 
CHO collaboration measurements by using our current 
clump sample in the fields analyzed previously. For 
example in the 24 fields used by Alcock et al. (1997a), 
we find 37 clump events giving r = 1.36+Q28 x 10 -6 , 
which is two sigma below the earlier result of r = 
3.9±i J x 10~ 6 . Over the 8 fields in our difference 
imaging study (Alcock et al. 2000a), our current, al- 
most independent sample gives t = 2.20lo 52 x 10~ 6 
at (I, b) = (2?68, — 3?35), which is in very good agree- 
ment with that earlier result of r = 2A3to^t x lcr<5 
(the agreement would be somewhat worse if we intro- 
duced the fudge factor discussed in the Introduction). 

Our new value of r = 2.17±£;^ x 10~ 6 is also 
very consistent with the recent EROS collaboration 
(Afonso et al. 2003) measurement. They found t = 
(0.94 ± 0.29) x 10~ 6 at (I, b) = (2?5, -4?0) from 16 



clump giant events and an efficiency calculation sim- 
ilar to ours. The values are different but they are re- 
ported at different position. The difference becomes 
entirely insignificant if we take into account either the 
numerical value of the gradient reported in §5 or read 
off the approximate optical depth value at b = —4.0 
from Figure 13. In addition, our new optical depth is in 
very good agreement with recent theoretical prediction 
of the Galactic bulge optical depth (t = 1.63 x 10~ 6 ) 
by Han & Gould (2003). 

We note that about 41% of the optical depth is in 
the events longer than 100 days (10 out of 62 events), 
and about 63% of the optical depth is in events longer 
than 50 days (21 out of 62 events). This may be at 
odds with some models of the Galactic structure and 
kinematics (e.g., Han & Gould 1996), but consistent 
with others (Evans & Belokurov 2002) which include 
streaming motion along the bar. The distributions of 
event durations are shown in Figure 14 where we plot 
histograms for uncorrected and efficiency-corrected 
cases. We note that the mean Einstein diameter cross- 
ing time of all our clump giant events is (i) — 56 ± 64 
days and for our fields close to the Galactic center 
(t) CGR = 39 ± 31 days. To facilitate proper com- 
parison with theoretical models we have to weigh each 
event by its inverse efficiency (Table 3). Then the av- 
erages become: (t) (eff) = 40 ± 50 days for all fields 
and (i) CGR (eff) = 30 ± 29 days for fields in the 
CGR 27 . Evans & Belokurov (2002) find average Ein- 
stein crossing times range from around 100 days in 
models that include bar streaming motion to 30 days 
in models without bar streaming. Our values are too 
uncertain to distinguish between these kinematic pos- 
sibilities and are consistent with both types of models. 

Note added. After this paper was completed, Sumi 
et al. (2005) presented determination of the microlens- 
ing optical depth toward the Galactic bulge based on 
OGLE-II data. Using a sample of 32 events from 
the extended red clump giant region of the color- 
magnitude diagram, they obtained r = 2.37tg 4g x 
10~ 6 at (l,b) = (l?16,-2?75). Their optical depth 
is fully consistent with our value of r = 2.17^0 38 x 
10~ 6 at (I, b) = (1?50, -2?68) which is based on 42 
clump giant events. 

We thank the referee, Przemyslaw Wozniak, for 



Because the distribution is non-Gaussian we also give the median 
and quartiles, which are (23.7, 12.0, 47.8) for all clump events and 
(18.7, 9.5, 46.9) for clump events in the CGR. 
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Fig. 6. — The upper panel shows the event with no 
color signal, and the lower panel the event with a weak 
signal. On each panel, the upper box shows the ampli- 
fication in the blue filter (the baseline flux has A = 1), 
the middle box the amplification in the red filter, and 
the lower box the difference between the two which is 
an indication of color. Shown are only the peak regions 
within t of the time with maximum amplification to- 
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Fig. 7. — Ratios of AV/AR for the 53 events with 
observations in the peak of the lightcurve. 
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Fig. 8. — Cumulative distribution of impact parame- 
ters weighted by inverse efficiencies (renormalized to 
55, the total number of unique non-binary clump giant 
events). Our cuts limit the impact parameter to 0.826. 
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Fig. 9. — Upper panel presents blue blend fraction dis- 
tribution for the clump sample from CGR (as defined 
in §5). The lower panel shows the recovered blend 
fractions for the sample of simulated unblended events 
(solid line) and events with input blend fraction of 0.5 
(dashed line). 
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Fig. 13. — Average optical depth in latitude and lon- 
gitude strips for events within 5.5° from the Galactic 
center. Gradient in b is much steeper than the one in I. 



Fig. 10. — Upper panel presents the distribution of S 
in blue filter for the clump sample from CGR (selected 
in the same way as for Figure 9). The lower panel 
shows the corresponding histogram for the simulated 
unblended events. 
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Fig. 11. — Clump detection efficiencies versus i for all fields, (a) solid green - 108, 113, 118, 1 19; red long dashed - 
1 16, 124, 125, 136, 142, 143, 146, 148, 149, 155; blue short dashed - 1 10, 159, 161, 162 (b) solid green - 1 1 1, 122, 
131, 132, 137, 150, 152, 158, 163, 171, 174; red dashed - 101, 104, 105, 109, 114, 120, 128 (c) solid green - 127, 
133, 138, 144, 147, 151, 153, 154, 156, 157, 160, 166, 168; red long dashed - 102, 103, 115, 121, 167; blue short 
dashed - 106, 107, 112, 117, 123, 126, 129, 130, 135, 141, 164, 165, 169, 170 (d) solid green - 176-180, 301-311, 
401-403; red dashed - 134, 139, 140, 145, 172, 173, 175. Note that fields 176-180 and fields 401-403 efficiencies are 
low because they were not observed for the whole 7 year period. As explained earlier, fields 301-3 1 1 are not included 
in the analysis but are included here for completeness. , „ 
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Fig. 12. — (a) Microlensing optical depth (xlO 6 ) toward each field that contains a clump event, (b) Optical depth in 
each field with a detected clump event and the number of events in that field. Only analyzed fields are shown in both 
panels. 




Fig. 14. — Left panel shows the histogram of the distribution of Einstein crossing times. The right panel shows the 
same distribution with each event weighted by its inverse efficiency and renormalized. 
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A. Theoretical Blending Estimates 

A.l. Fraction of Faint Star Blend Events 

A simple theoretical estimate of the effect of blending on our clump giant analysis can be made with the help of 
the HST luminosity function for Baade's Window from Holtzman et al. (1998). We multiply this luminosity function 
by a factor of 1 .5 to account for the fact that our Central Galactic Region (CGR; defined in §5) is at a lower Galactic 
latitude than Baade's Window. Faint stars within a radius of approximately full width at half maximum of the point 
spread function of an image will be close enough to the clump giant star to be blended with it. Both the baseline 
level and the microlensing signal for each event are dominated by images with seeing better than the median, so we 
make the conservative choice of the median seeing in the MACHO blue band images, namely 2.1", to use for our 
blending estimates. The Holtzman et al. (1998) luminosity function multiplied by 1.5 implies a total of <~ 16 stars 
with an unreddened absolute V magnitude My > 1.3 within our seeing disk, and combined brightness of these stars 
is My = 3.13, or 2.13 magnitudes fainter than a typical clump giant. Due to turnoff feature, My = 3.13 is too bright 
for a single bulge main sequence star, but three My — 4.32 stars could combine to give the same total brightness. 
The blending effect of these stars will be approximately equivalent to the blending effect from the luminosity function 
under consideration. 

We consider an object that consists of a single clump giant with blend fraction f cg = 0.877 blended with three 
main sequence stars, each with blend fraction, / ms = 0.041. In order to pass our detection threshold of A^ x > 1.5 
for fits without blending (apparent observed) which corresponds to u°^t = 0.827, the blended clump giant must reach 
an actual brightness of ^4 max , cg > 1.57, corresponding ?i cu t, C g = 0.771. For each of the three faint, main sequence 
stars, the requirement of A^ x > 1.5 implies a peak magnification of A maK ms = 13.2 or u cu t. m s = 0.076. The 
blending will also cause a reduction in the fit t values for a fit without blending compared to the actual values. While 
the threshold magnification is > 1.5 for detection, the typical event has A^l x — 2.5, so we will use this 

magnification for our estimates of the blending effect on t. For the clump giant, the unblended fit will yield a t° g s 
value about a factor of 1 .06 smaller than actual value, t cg and for the faint main sequence stars, the true t ms value will 
be reduced by a factor of 8.2 according to a calculation similar to that shown in Fig. 2 of Wozniak & Paczyhski (1997). 
This results in a negligible change in the detection efficiency for the clump giant by a factor of £ cg = 1.00, but for a 
main sequence star event with a true t ms <~ 50 days, the detection efficiency is reduced by about a factor of £ ms = 1.5. 
Using the fact that microlensing events are uniformly distributed in u m ; n , we can find the ratio of the numbers of faint 
blended events to clump giant events, 

^ms _ 3 Mcut.ms/ Cms _ q (Al) 
-^cg ^cut,cg/ £cg 

So, this estimate indicates that <~ 83.5% of apparent clump giant events should have actual clump giant sources, while 
<~ 16.5 % of the apparent clump giant events will be due to much fainter main sequence stars. This estimate has not 
included a contribution from lensing of faint blended stars that are part of an apparent clump giant that is actually a 
blend of stars that are all fainter than a clump giant. Such events are rare because of the steep drop in the luminosity 
function just below the clump and the fact that most of the stars that are slightly fainter than the clump are much bluer 
foreground disk main sequence stars that are much less likely to be lensed. So, it is unlikely that the fraction of blended 
events could reach 25%. 

A.2. Blending Effect on Microlensing Optical Depth 

The general formula for the microlensing optical depth is 

Advents f 
_ 7T ^ ti 



(A2) 



where E = iVT Q b s is the total exposure which depends on the number of source stars, N, and the time span of the 
observations, T b s - Let us assume that the number of clump giants in the bulge fields have been estimated properly 
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and so N does not have to be analyzed in what follows, and T Q b s is fixed by the experimental setup. Now, we seek to 
calculate the correct optical depth for non-contaminated clump giant sample that would have been observed if there 
had been no blending. 



T *rue 



^would-be 



CS 4£ g e(C)' (A3) 

where N™° uld ~ hc is the number of clump events that would have been observed if there had been no blending. Note 
that: 

Arwould-bc obs 

_cg _ _^cut (A4) 



-^cg ^cut,cg/£cg 

where typically we do not expect AT cg = N cvcnts (see equations Al and A6). 

In reality, we effectively evaluate the following estimator in the clump region of the color-magnitude diagram: 

^cg Xobs ^ms Xobs 

obs = obs , obs = Jl_ c S,i i 71 V 'msj , ^ 

where r°g s and r^ s are the contributions to the observed optical depth coming from true lightly-blended clump giants 
and from main sequence blends, respectively. Therefore: 

Ncg + N ms - Advents, (A6) 

and we know from equation (Al) that: 



au = - at cs . (a?) 



3 w cu t ims / 

1*cut,cg/ £cg 

Let us introduce the following notation: r] cg = i° g s /i C g and ?7 ms = i^/i ms . Then we have: 



true _ n U cut AT I teg \ (AQ\ 

obs ^ £cg I t cg \ 

„obs _ ^ ^ms 3 Mcut,ms/£ms Ar / ''ma , (A 10) 



ms 4E Vms Ucut ,cg/ec g S \e(£ m s)/' 
where notation (.) indicates an average over events. To make further progress we will assume that: 



tcs (All) 



(^ms) / \ e (^cg) 

but the conclusions will be almost identical even if equation (All) does not hold precisely. Combining equations (A8), 
( A9), (A 1 0), ( A5), and (A 1 1 ) we find: 

1 3 

b s ^cut,cg ^cut.ms 

Z = Z^s ^ (A12) 

,-true 7/ obs ^ ' 

eg u cut 

We can apply equation (A 12) to the situation from § A.l, where we considered the effect of blending on typical 
events with A^ x w 2.5. For the estimate of the fraction of blended events, we implicitly assumed that the blended 
clump giant and main sequence stars all stand out above the background, so that the photometry code identifies a 
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blended star that has the combined brightness of the clump giant and the main sequence stars. However, the unresolved 
main sequence stars also contribute to the background brightness that is attributed to the sky background by the 
photometry code. In areas of the bulge fields where giant stars are absent, our photometry code finds "objects" that are 
generally brighter than the top of the bulge main sequence, and blends of multiple main sequence stars. Thus, the faint 
main sequence stars contribute to the "sky background" if they occur in regions of lower than average star density, and 
they can contribute to stellar "objects" identified by our photometry code if they are in high density regions. When a 
clump giant is present, it will be blended with both the faint stars that make up part of the background and with the 
stars that would have combined to make up part of a stellar "object" if the clump giant happens to be in a location with 
a higher than average main sequence star density. So, to bracket the real situation, we will consider two cases in which 
the clump giant is blended with the same 3 main sequence stars discussed in § A. 1 : 

1. The faint main sequence stars and the clump giant will combine to form a single stellar "object" with their 
combined brightness. 

2. The faint main sequence stars will only contribute to the background "sky" brightness, and the stellar "object" 
seen by the photometry code will have the same brightness as the clump giant. 

In case 1, i reduction factors for the fits assuming no blending are r] cg = 1.06 and ry ms = 8.2, whereas the u min 
thresholds are u cu t,cg = 0.771, u cu t,ms = 0.076, and !i* t s = 0.827. Therefore, using equation (A12) we obtain: 



T obs 



= 0.913. (A13) 



In case 2, t reduction factors for the fits assuming no blending are rj cg = 1.0 and ?7 ms = 7.4, whereas the u m i n 
thresholds are w cut , C g = «°ut = 0-827 and u C ut,ms = 0.086. Therefore, equation (A12) yields: 

_obs 

1.042. (A14) 



—true 

eg 



Thus, our two estimates of r obs bracket the true value. 

There are some additional blending effects that we did not include in this estimate. There will be some "objects" 
that pass our clump giant cuts only because of blending. Typically, these will be dominated by light from a single star 
that is slightly below the clump giant cut, and the inclusion of these false clump giants will be partly compensated by 
the loss of true clump giants that appear brighter than our clump region due to blending. Overall, we expect this effect 
to cause a slight increase in T obs /r*g UC , but the overall systematic error in the optical depth due to blending is probably 
less than 10%. 

The selection of the clump giant sample minimizes the effects of blending on the optical depth measurement because 
clump giants are primarily blended with stars that are much fainter. This means that blending has only a small effect 
on the apparent t° g s values for the true clump giant events that dominate the sample, while the faint star blends have 
such a large duration shift that they cause little bias of the microlensing optical depth. 



B. Optical depth estimation 

The optical depth is usually defined as the probability of a microlensing event happening at a given time on any 
given star. An "event" is usually defined as occurring when a lens comes within one Einstein radius of the line-of-sight 
to the source. This is easily determined observationally since the formula for magnification is 

M u ) = t4t^ > (b1) 

where u = I /Re is the distance of the lens from the line-of-sight in units of the Einstein radius Re, and 

R^2{ GMD ^- D ^ \' 2 . (B2, 
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where M is the mass of the lens, D s is the observer-source distance, and Di is the observer-lens distance. Thus u < 1 
or A > 1.34 means the line-of-sight is within one Einstein radius. 

To estimate the probability of lensing occurring we must use the observed microlensing events. First consider the 
case of a 100% efficient microlensing experiment. When there are many stars observed, one would like to find the 
probability of lensing per star at a given time by dividing the number of stars being lensed at a given time by the total 
number of observed stars. Observationally, however, one cannot tell whether or not a given star is being lensed; one 
only measures stars brightening and then fading as the line-of-sight enters and leaves the Einstein radius. In this case 
one can estimate the lensing probability by the fraction of the total available observing time that a given star spends 
lensed with A > 1.34. With many stars, the optical depth is thus estimated as r = ^iU/E, where ti is the time 
inside the Einstein radius (i.e. A > 1.34), and the exposure E = -ZV*T b s is the number of observed stars times the 
duration of the observations. The quantity ti is easily found from the two microlensing fit parameters: ti, the Einstein 
diameter crossing time, and A maXji , the maximum magnification: 

t i =i i ^/l-u 2 min ., (B3) 

where w m i n! i is found from A maXt i by inverting equation (Bl) above. 

For an imperfect observational program only a fraction e of microlensing events will be detected. Calculating this 
efficiency as a function of event duration and magnification allows a correction to be made for each detected event 
giving the formula we use: 

N events i 

If one wished to use the fit U and avoid using U one could replace ti in the formula above with an estimate of the 
typical value of t. This could be done by noting that every impact parameter, u, between and 1 is equally likely 
(given perfect efficiency) and so averaging over all impact parameters for a given Einstein ring gives (ti) = j (ti). 
Since for a large number of events the sum over ti is effectively an average over t i7 one could replace ti with jti inside 
the sum. The resulting formula is equation (12), and has been used in the past. 

One difference between the two formulas is that the new formula allows weighting each event with an efficiency 
calculated as a function of both ti and w m i n ,i- When using the sum over ti one usually weights by the efficiency only 
as a function of t. 

There is one subtlety when calculating efficiencies for each event, rather than as a function of t. Since we used a 
cut A max > 1.5 and optical depth is defined as the probability of finding events with A max > 1.34, the above formula 
needs to be multiplied by a correction factor of 1.09 to account for events with 1.34 < A max < 1.5. These low 
amplification events are not detected so can not be included in the efficiency calculation which only inputs simulated 
events for the measured values of t and A max . When efficiencies are calculated as a function of t, events with such 
values are input, but not detected, resulting in a lower efficiency (and therefore higher optical depth) for that value of 
t. One can invert equation (Bl) and find that u(1.5) = 0.83, and u(1.34) = 1. With the assumption that the lens is 
moving at uniform velocity through the line of sight, the correction factor becomes simply the ratio of area within the 
Einstein ring to the area in the ring with an impact parameter less than 0.83. In Figure 15 one can see the Einstein ring 
and the ring at u — 0.83 ( A max = 1.5). The areas missed with our cuts are the small sections on the top and bottom 
of the circle. 



C. Upper limits on optical depth 

Optical depth in a given field can be estimated as a number of events in a given field multiplied by the contribution 
of a typical event. Therefore, one can define upper limits on the optical depth in fields with no events as: 

lim.CL ,,CL 1 event <ni\ 
T f =V T f , (CI) 
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Fig. 15. — Area of Einstein rings used for calculating optical depths, shaded regions are missed by our cuts. 

where rj cvcnt would be a typical contribution of one event to the optical depth in field / and i/ CL is a multiplication 
factor dependent on the confidence level (CL). Equation (CI) can be expressed as: 

lim.CL CL n I 1 \ 

Tf =v ^f^Xwi (C2) 

where the angle brackets indicate an average in field /. Let us first note that for each field, we have efficiencies e(t), 
which come from our simulations. However, since we have no events in the fields of interest we cannot estimate the 
right side of equation (C2) based exclusively on the data in field /. Therefore, we are forced to assume that duration 
distribution in the entire area of our bulge survey is universal and proceed as follows. Based on the observed duration 
distribution in fields with events and efficiencies in those fields we can recover the intrinsic duration distribution of 
events. This distribution can then be mapped using the efficiencies in field / to recover the hypothetical duration dis- 
tribution that would be observed in field /. Each ti from our clump sample contributes to the hypothetical distribution 
in field / with a weight Wfj = €f(ti)/ei(ii), where e/ is the efficiency in field /, and ei is the efficiency in the event's 
field of origin. We now may express: 



t \ £f=7 nt8 ™/,; 



w / / Z^=i w f,i 



^2i=i U/ti(h) 
EilT ntS e/(*i)Ai(*0 
Substitution of equation (C3) to equation (C2) leads to: 

lim.CL _ CL n 12^=1^" ^i/^Si) ,^y^ 

which is the formula that we use in §5. Let us explain the role of the multiplication factor v clj . The number of 
observed events is Poisson distributed. To constrain optical depth we simply ask what is the upper limit on the number 
of expected events given detected events for the Poisson distribution. We can obtain the necessary number of 
expected events v ch that rules out T j lm ' CL or larger at CL from the following equation: 

, PTn 1.0 -CL 
exp(^ CL ) = (C5) 



For la (68.17%) CL, one obtains v lrT = 1.8379. 
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Table 1 
Selection Criteria 



Selection 



Description 



microlensing fit parameter cuts: 

N {V -r) > 

Xlut < 3.0, xLt > 0.0 



Namp > 8 
N rising 2^ 1 1 N 



falling ^ 1 
A m ax > 1-5 
(A max 

- 1) > 2.0((7R + (7b) 
(CTfl + ffs) < (0.05 5 X 2 )/(^Vam P 

(Namp xLt fchrom)/(befaftS X 2 ) < 0.0003 
(N amp xl n )/{befaftS X 2 ) < 0.00004 
5 X 2 /fc2 > 320.0 

0.5(rcrda + bcrda) < 143.0 



require color information 

require high quality baselines 

require 8 points in the amplified region 

require at least one point in the rising and falling part of the peak 

magnification threshold 

signal to noise cut on amplification 

remove spurious photometric signals caused by nearby saturated stars 

good overall fit to microlensing lc shape 
same quality fit in peak as whole lc 
source star not too crowded 



£& ut 7£&"*° < 2.0 

t > 419.0, t < 2850.0 

i < 1700 


remove long period variables 

constrain the peak to period of observations 

limit event duration to ^half the span of observations 


clump giant cuts: 




V > 15, V < 20.5 


select bright stars with reliabale photometry 


V > 4.2(V - R) + 12.4 


define bright boundary of extinction strip 


V < 4.2(V - R) + 14.2 


define faint boundry of extinction strip 


(V — R) > (V — R)boundary 


avoid main sequence contamination 


exclude fields 301 — 311 


avoid disk contamination 



NOTE. — Designations V and (V — R) indicate baseline quantities, i.e. the ones in the limit of no microlensing-induced 
amplification. 
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Table 2 
Event parameters 



field.tile.seq 


RA (J2000) 


DEC (J2000) 


I 


b 


V 


(V — R) 


(V R~) boundary 


to 


t 


j4max 


X 
n dof 


101.21689.315 


18 06 58.32 


-27 27 45.7 


3.639 


-3.327 


17.50 


1.09 


0.81 


593.5 


165 


4.80 


1.19 


102.22466.140 


18 08 47.04 


-27 40 47.3 


3.643 


-3.783 


16.81 


0.95 


0.81 


1275.6 


153 


3.88 


0.95 


104.19992.858 


18 02 53.76 


-27 57 50.0 


2.761 


-2.784 


18.59 


1.09 


0.85 


1169.7 


106 


2.23 


2.15 


104.20119.6312" 


18 03 20.64 


-28 08 14.6 


2.658 


-2.955 


17.99 


1.06 


0.78 


1975.3 


38.7 


1.62 


1.62 


104.20251.50 


18 03 34.08 


-28 00 19.1 


2.797 


-2.933 


17.33 


0.93 


0.82 


493.9 


287 


10.1 


1.75 


104.20251.1117 


18 03 29.04 


-28 00 31.0 


2.785 


-2.919 


18.22 


0.98 


0.82 


463.2 


45.5 


2.09 


0.93 


104.20259.572 


18 03 33.36 


-27 27 47.2 


3.269 


-2.665 


18.16 


1.19 


0.91 


538.5 


14.2 


1.58 


1.56 


104.20382.803 


18 03 53.28 


-27 57 35.6 


2.872 


-2.973 


17.79 


0.96 


0.84 


1765.8 


254 


5.72 


4.40 


104.20515.498 


18 04 09.60 


-27 44 35.2 


3.090 


-2.919 


17.63 


0.99 


0.81 


2061.7 


53.6 


1.76 


1.42 


104.20640.8423 6 


18 04 33.60 


-28 07 31.8 


2.800 


-3.183 


17.15 


0.98 


0.79 


2374.9 


30.8 


1.65 


2.23 


104.20645.3129 


18 04 26.16 


-27 47 35.2 


3.077 


-2.997 


17.68 


1.07 


0.83 


1558.3 


152 


1.82 


0.77 


105.21813.2516 


18 07 06.96 


-27 52 34.3 


3.292 


-3.555 


17.41 


1.08 


0.83 


1928.5 


17.5 


13.2 


1.14 


108.18947.3618 


18 00 29.04 


-28 19 20.3 


2.186 


-2.498 


17.72 


1.10 


0.89 


1287.9 


45.6 


2.05 


0.94 


108.18951.593* 


18 00 33.84 


-28 01 10.6 


2.458 


-2.363 


18.07 


1.18 


0.95 


1991.6 


47.3 


4.05 


7.02 


108.18951.1221 


18 00 25.92 


-28 02 35.2 


2.423 


-2.350 


18.94 


1.18 


0.95 


583.6 


45.8 


2.16 


0.92 


108.18952.941 


18 00 36.00 


-27 58 30.0 


2.501 


-2.348 


18.90 


1.16 


1.02 


1324.9 


66.9 


2.39 


1.11 


108.19074.550 


18 00 52.32 


-28 29 52.1 


2.076 


-2.659 


17.54 


1.03 


0.88 


2050.6 


9.30 


1.78 


0.97 


108.19334.1583 


18 01 26.40 


-28 31 14.2 


2.117 


-2.779 


17.67 


0.94 


0.86 


1231.4 


18.5 


4.09 


0.38 


109.20119.1051 a 


18 03 20.64 


-28 08 14.6 


2.658 


-2.955 


18.09 


0.97 


0.78 


1977.6 


28.7 


1.62 


1.01 


109.20640.360 6 


18 04 33.60 


-28 07 32.5 


2.799 


-3.183 


17.33 


1.05 


0.79 


2375.1 


28.7 


1.67 


0.72 


110.22455.842 


18 08 51.36 


-28 27 11.2 


2.971 


-4.169 


18.28 


0.98 


0.76 


1900.1 


19.8 


2.54 


0.84 


113.18552.581 


17 59 36.00 


-28 36 24.1 


1.842 


-2.470 


17.60 


1.02 


0.86 


603.7 


27.9 


1.77 


1.11 


113.18804.1061 


18 00 03.36 


-29 11 04.2 


1.390 


-2.843 


18.28 


1.00 


0.84 


1166.6 


6.50 


1.62 


0.96 


113.19192.365 


18 01 06.96 


-29 18 53.6 


1.391 


-3.109 


17.80 


1.17 


0.96 


2092.6 


36.9 


1.97 


0.52 


114.19712.813 


18 02 11.52 


-29 19 20.6 


1.500 


-3.317 


18.12 


1.08 


0.87 


1643.4 


23.7 


3.58 


0.62 


114.19846.777 


18 02 36.72 


-29 01 41.9 


1.801 


-3.252 


17.98 


1.02 


0.80 


545.4 


65.0 


1.55 


1.86 


114.19970.843 


18 02 54.48 


-29 26 29.4 


1.472 


-3.511 


17.99 


0.96 


0.85 


1281.9 


14.7 


4.47 


0.65 


118.18014.320 


17 58 25.20 


-29 47 59.6 


0.678 


-2.840 


17.01 


0.97 


0.89 


877.6 


25.5 


1.58 


1.05 


118.18141.731* 


17 58 36.72 


-30 02 19.3 


0.491 


-2.995 


17.92 


1.06 


0.94 


884.1 


9.80 


12.6 


0.65 


118.18271.738 


17 58 52.56 


-30 02 08.2 


0.522 


-3.043 


18.25 


1.18 


0.96 


456.4 


53.7 


2.23 


0.69 


118.18402.495 


17 59 13.92 


-29 55 53.4 


0.651 


-3.058 


17.63 


1.08 


0.93 


454.1 


23.8 


3.60 


0.96 


118.18797.1397 


18 00 06.96 


-29 38 06.0 


1.004 


-3.078 


19.34 


1.30 


0.97 


2013.1 


126 


8.35 


1.07 


118.19182.891 


18 01 09.84 


-29 56 19.0 


0.852 


-3.425 


17.97 


0.93 


0.84 


2367.0 


13.6 


6.28 


1.20 


118.19184.939 


18 01 10.32 


-29 48 55.4 


0.960 


-3.366 


18.07 


0.99 


0.83 


2692.8 


74.4 


2.14 


1.13 


121.22032.133 


18 07 46.80 


-30 39 41.4 


0.915 


-5.025 


16.37 


0.86 


0.75 


1287.3 


19.5 


1.58 


2.51 



TABLE 2 — Continued 



field.tile.seq 


RA (J2000) 


DEC (J2000) 


I 


b 


V 


(V-R) 


(V ^)boundary 


to 


i 


^max 


2 

X 

aoj 


158.27444.129 


18 20 10.32 


—25 08 25.8 


7.097 


-4.827 


16.97 


0.76 


0.63 


2376.4 


41.7 


2.31 


1.80 


162.25865.442 


18 16 30.72 


—26 26 58.2 


5.549 


—4.712 


17.52 


0.90 


0.73 


1225.3 


60.9 


1.58 


1.32 


162.25868.405 


18 16 46.08 


-26 11 43.1 


5.801 


-4.643 


17.33 


0.83 


0.77 


1222.8 


30.8 


2.77 


1.39 


176.18826.909 


18 00 15.36 


—27 42 47.9 


2.691 


-2.153 


18.77 


1.20 


1.01 


2364.7 


24.3 


1.77 


1.17 


178.23531.931 


18 11 17.52 


—25 59 51.4 


5.391 


—3.467 


18.71 


1.10 


0.92 


1590.3 


13.7 


2.18 


1.08 


180.22240.202 


18 07 57.84 


— 25 24 48.2 


5.542 


—2.528 


18.92 


1.40 


0.98 


1286.9 


311 


3.42 


1.15 


401.47991.1840 


17 57 07.44 


— 28 13 44.0 


1.899 


— 1.810 


19.46 


1.37 


1.12 


1953.0 


44.7 


1.92 


0.61 


401.47994.1182 


17 57 07.92 


f\Ci f\f\ ^ "~T "~T 

— 28 00 27.7 


2.091 


—1.701 


19.65 


1.53 


1.23 


2362.7 


60.6 


3.68 


1.05 


401.48052.861 


17 57 23.76 


— 28 10 24.6 


1.977 


— 1.834 


18.75 


1.31 


1.13 


2291.5 


72.0 


2.65 


0.69 


401. 48167.1934 


17 57 57.36 


— 28 27 09.7 


1.796 


—2.080 


18.42 


1.23 


1.09 


2179.0 


125 


2.53 


0.89 


401. 48229. 760 + 


17 58 13.44 


—28 20 32.3 


1.921 


—2.076 


18.07 


1.15 


1.00 


2694.5 


20.1 


1.82 


1.18 


401.48408.649* 


17 59 08.88 


-28 24 54.7 


1.959 


—2.289 


17.70 


1.20 


0.92 


2325.4 


74.6 


3.07 


1.30 


401.48469.789 


17 59 22.32 


-28 20 21.8 


2.050 


—2.294 


18.23 


1.17 


0.92 


1554.8 


5.70 


5.96 


0.42 


402.47678.1666 


17 55 22.56 


— 29 04 12.7 


0.978 


— 1.901 


19.19 


1.31 


1.07 


2393.1 


13.6 


2.48 


0.88 


402.47737.1590 


17 55 41.52 


-29 09 00.7 


0.944 


—2.001 


18.49 


1.04 


0.99 


1552.3 


5.50 


2.79 


1.03 


402.47742.3318 


17 55 42.96 


—28 49 17.4 


1.230 


— 1.840 


20.19 


1.62 


1.10 


1377.1 


50.1 


8.64 


1.21 


402.47796.1893 


17 56 10.80 


rtA 1 i~\ A A A 

— 29 10 44.4 


0.972 


—2.107 


19.52 


1.33 


1.03 


1678.7 


19.3 


22.5 


1.47 


402.47798.1259 


17 56 11.04 


—29 06 14.0 


1.038 


-2.070 


19.72 


1.36 


1.10 


1583.7 


47.6 


4.94 


1.38 


402.47799.1736 


17 56 07.20 


-28 58 32.5 


1.142 


— 1.994 


19.59 


1.29 


1.16 


2410.0 


115 


2.26 


0.47 


402.47856.561 


17 56 24.96 


-29 12 56.5 


0.966 


-2.170 


18.81 


1.38 


1.15 


1566.6 


12.3 


2.14 


1.63 


402.47862.1576* 


17 56 20.64 


-28 47 42.0 


1.323 


-1.945 


19.47 


1.30 


1.15 


2028.5 


54.3 


7.52 


6.13 


402.48158.1296 


17 57 47.52 


-29 03 54.4 


1.247 


-2.346 


18.69 


1.12 


1.01 


2325.7 


15.9 


3.94 


0.98 


402.48280.502 


17 58 24.96 


-28 57 46.4 


1.404 


-2.422 


18.17 


1.18 


0.90 


1342.5 


25.5 


2.08 


1.14 


403.47491.770 


17 54 38.64 


-29 33 13.0 


0.480 


-2.007 


17.94 


1.29 


1.04 


2633.9 


48.5 


2.43 


0.99 


403.47550.807 


17 55 00.00 


-29 35 03.8 


0.492 


-2.089 


17.95 


1.21 


1.04 


2774.3 


11.6 


3.60 


0.87 


403.47610.576 


17 55 17.04 


-29 37 40.8 


0.486 


-2.164 


17.54 


1.13 


0.98 


2657.9 


9.60 


4.75 


2.56 


403.47845.495 


17 56 22.80 


-29 55 16.3 


0.351 


-2.517 


18.20 


1.24 


1.07 


1945.8 


15.0 


2.03 


0.88 



NOTE. — Designations: 
* : binary event, 

a and b : each superscript marks two members of a pair that represents a single microlensing event in two overlapping fields. 



Table 3 

Optical Depth Contribution by Event 



field. tile. seq 


l i 


i 


A 




e(t) 




ti) 




101. 


.21689 


.315 


161 


165 


4.80 


0. 


.68 


0.55 


1 


.62 


1.47 


102. 


.22466 


.140 


147 


153 


3.88 





.39 


0.32 


3 


.67 


3.34 


104. 


1 9999 


858 

-OfJO 


93.4 


107 


2.23 


n 


55 


0.43 


1 


.12 


1.19 


104. 


.20119 


631 2 


25.7 


38.3 


1.62 


n 


48 

■ 4:0 


0.38 





.36 


0.47 


104. 


.20251 


50 


285 


287 


10.1 


n 

VJ . 


77 

. I I 


0.63 


2 


.44 


2.18 


104. 


.20251 


.ILL i 


38.7 


45.6 


2.09 


VJ. 


46 

4:VJ 


0.39 





.56 


0.56 


104 


20259 


579 

• O 1 


9.22 


14.2 


1.58 


(1 

U ■ 


.41 


0.33 





.15 


0.20 


104 


90389 


803 

.OVJO 


250 


254 


5.73 


n 

VJ . 


78 

■ 1 O 


0.62 


2 


.11 


1.95 


104. 


.20515 


.498 


40.4 


53.5 


1.76 


0. 


.48 


0.39 





.56 


0.65 


104. 


.20640 


.8423 


21.4 


30.8 


1.65 


0. 


.47 


0.4 





.30 


0.37 


104. 


.20645 


.3129 


118 


152 


1.82 


0. 


.66 


0.52 


1 


.19 


1.38 


105. 


.21813 


.2516 


17.5 


17.5 


13.2 


0. 


.40 


0.31 





.35 


0.32 


1 08 

-LVJO . 


1 8947 


361 8 


38.4 


45.6 


2.05 


n 


56 


0.48 





.38 


0.38 


1 08 

-LVJO 


1 8951 

. 1-OUO -L 


593 


45.8 


47.3 


4.05 


n 

VJ. 


57 


0.47 





.44 


0.39 


1 08 

-LVJO ■ 


1 8951 

. -L OUO J- 


.1221 


39.5 


45.8 


2.16 


n 


57 


0.48 





.38 


0.38 


-LVJO ■ 




Q41 


59.7 


66.9 


2.39 


n 

VJ . 


■ oo 


0.46 





.56 


0.57 


1 OR 

-LVJO . 


1 Q074 


550 

■ OOO 


7.07 


9.26 


1.78 


n 

VJ . 


.41 


0.34 





.09 


0.11 


1 08 

-LVJO ■ 


1 9334 


1 583 

. -LfJOO 


17.9 


18.5 


4.09 


n 


53 


0.42 





.18 


0.17 


1 OQ 

-LVJi? ■ 


.20119 


1 051 

. -LVJtJ -L 


19.4 


28.7 


1.62 


n 

VJ . 


.44 


0.4 





.29 


0.34 


109. 


.20640 


.360 


20.4 


28.7 


1.67 


0. 


.45 


0.4 





.29 


0.34 


110. 


.22455 


.842 


18.0 


19.8 


2.54 


0. 


.41 


0.33 





.46 


0.46 


113 


.18552 


.581 


21.2 


27.9 


1.77 


0. 


,53 


0.46 





.21 


0.23 


1 1 3 

_L lu . 


1 8804 

. J_ OOVJ4:. 


1 061 

- -LVJU -L 


4.44 


6.53 


1.62 


n 

vj. 


29 


0.27 





.08 


0.09 


1 1 3 

_L lO . 


.19192 


365 

• OOO 


30.4 


36.9 


1.97 


n 


59 

■ OU 


0.49 





.27 


0.28 


114 


1 971 2 


81 3 


22.7 


23.7 


3.57 


n 

u . 


.44 


0.36 





.36 


0.33 


114. 


1 9846 

. _L C/04:VJ 


777 

• III 


40.0 


64.9 


1.55 


n 


46 


0.39 





.61 


0.84 


114 


1 9970 

. -L U U t VJ 


843 

.04:0 


14.3 


14.7 


4.47 


n 

VJ . 


45 
^±o 


0.32 





.22 


0.23 


118. 


.18014 


.320 


16.3 


25.5 


1.58 


0. 


.54 


0.44 





.18 


0.24 


118. 


.18141. 


.731 


9.74 


9.77 


12.6 


0. 


.45 


0.34 





.12 


0.12 


118 


.18271 


.738 


46.9 


53.7 


2.23 





.56 


0.46 





.49 


0.49 


118 
±±o. 


1 8409 




22.8 


23.8 


3.60 


n 
u. 


^4 


0.43 





.24 


0.23 


118. 


.18797 


.1397 


125 


126 


8.35 





.73 


0.58 


1 


.00 


0.91 


118. 


.19182 


.891 


13.4 


13.6 


6.28 


0. 


.49 


0.39 





.16 


0.15 
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0. 
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TABLE 3 — Continued 



field. tile. seq 


k 


t 


A 




e(t) 


r(U) 


r(t) 
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.48229 
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20.1 


1.82 


0.24 
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0.26 
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.807 


11.1 


11.6 


3.60 


0.21 


0.16 


0.20 


0.20 


403 


.47610 


.576 


9.37 


9.59 


4.76 


0.2 


0.14 


0.19 


0.19 


403 


.47845 


.495 


12.6 


15.0 


2.03 


0.25 


0.20 


0.20 


0.22 
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Table 4 

Alternative designations of our events 



field.tile.seq RA (J2000) DEC (J2000) MACHO Alert 



EROS 



OGLE 



i m 
1U1 


Ol RQCl 


Q1 K 

olo 


1 Q 

lo 


UO 


Jft 


"30 
.JZ 


117 

— Z / 


0~7 

Z / 


A C 

4j 


n 
. 1 


1 n 1 d 

1UJ-D 




i no 
lUz 


OO/I (Zfc 
.ZZ400 


1 AP 
14U 


1 
1 ft 


no 
Uo 


4/ 


.U4 


07 


/in 
4U 


An 

4 /. 


-i 

.3 


QC DT 

yj-rSlAj- 


1 ^ 
- 1 j 


1 C\A 

1U4 


i nnno 


oro 

ooo 


1 
1 ft 


no 
Uz 


J J 


"7A 


— Z ( 


j / 


cn 

jU. 


n 
U 






1 C\A 

1U4 


om i n 

.zuiiy. 


OOlZ 


1 
1 ft 


U3 


on 
ZU 


A/1 
.04 


no 
— 28 


nQ 
Uo 


1 A 

14 


.0 


m DT r 


1.A 

■34 


1 PA 


ono c: i 
.zUzol 


en 


1 O 


U3 


34 


no 
-Uo 


no 
— 25 


nn 
UU 


i n 

iy 


1 

, J 


1 n/i 
1U4-L. 




1 PA 
1U4 


ono^i 
-ZUZOl 


1117 

111 / 


1 
1 ft 


n^ 
U3 


OO 

zv 


n/1 
.U4 


OS 

— 2o 


nn 
UU 


J J. 


n 

,U 


1U4-JJ 




1 PA 

1U4 


-zUzoy. 


( z 


1 
1 ft 


no. 
U3 


j j 


Q A 
.JO 


— 2 / 


07 
z / 


An 
4/. 


.z 


1 O/l A 

1U4-A 




1U4 


onQ QO 
.zUooz. 


ono 
oUo 


1 O 
1 ft 


n^ 
U3 


J J 


O Q 
.Zft 


n7 

— 2 / 


^7 

J / 


jj. 


,o 


nA DT Z" 1 


1 o 
■Iz 


1 C\A 

1U4 


on cci c 
.zUolo. 


a no 
4yo 


1 

1 ft 


n,i 
U4 


no 

uv 


An 
,0U 


n7 

— 2 / 


A A 

44 


jj. 


.z 


m DT r 


C 

-Jft 


i n a 
1U4 


one a n 


O A OO 

,o4zo 


1 o 

1ft 


f\A 

U4 


JJ 


An 
.oU 


OO 

— 20 


U/ 


31 


o 
.ft 






i n a 
1U4 


one? a c: 


o i on 

oizy 


1 o 

lo 


ni 
U4 


ZD 


1 A 
.10 


— 2 / 


,-1 "7 

4/ 


JJ. 


,z 


nA DT z -1 


1 

- 1 


1 PK 
1U0 


Ol Q1 Q 

-zlolo. 


ZOlO 


1 
1 ft 


U / 


UO 


OA 

.yo 


OT 7 

— 2 / 


JZ 


'J/l 
J4 


.3 


n7 dt r 
y /-rSlAj- 


i n 
- JU 


i no 
lUo 


.ioy4 / . 


oolo 


1 O 
1 ft 


nn 


on 
ZV 


Pi A 

.U4 


OO 

— 28 


i n 
IV 


on 

zU. 


-3 
. J 






i no 
lUo 


1 one i 

.loyoi. 


cnQ 

oyo 


1 

1 ft 


nn 
UU 


j J 


Q A 

,o4 


OO 

— 28 


n 1 
UJ 


1 n 

LU. 


.0 


m dt r 


OO 

■Zo 


i no 
lUo 


i one: i 

.loyoi 


1 OOI 

Izzl 


1 

1 ft 


nn 
UU 


O 5 

ZJ 


no 
.Vz 


OO 

— 28 


no 
Uz 


jj. 


.z 


1 no a 
lUo-A 




i no 
lUo 


.loyoz. 


n /i i 

y4i 


1 o 

1ft 


nn 
UU 


JO 


nn 
.UU 


07 

— 2 1 


CO 

Jft 


Qn 
JU. 


.U 


nc DT 


i j 

■JZ 


i no 
lUo 


1 PPT A 

.iyu ( 4 


DOU 


1 
1 ft 


nn 
UU 


JZ. 


QO 
. JZ 


OO 

— 28 


OO 

zv 


CO 

JZ. 


, J 


n7 dt r 
y /-dIAj- 


CO 

oy 


i no 
lUo 


1 OQQ A 

.iyoo4. 


1 KOQ 

looo 


1 

1 ft 


n i 
UJ 


OA 
ZD 


/i n 
.4U 


OO 

— 28 


J 1 


1 A 

14. 


-> 
.z 


QC dt 


1 A 

-14 


i nn 

iuy 


om i n 
.zUliy 


lUol 


1 

1 ft 


n^i 
Uj 


on 
zU 


A/l 

.04 


OO 

— 28 


no 
Uo 


1 /i 
14. 


O 






iuy 


.ZU04U 


Q^n 
ooU 


1 
1 ft 


n/i 
U4 


j j 


An 
.oU 


OS 

— 28 


07 

u / 


70 

JZ. 


c 

. J 






i 1 n 
11U 


OO A £ K 

.zz4oo 


O A O 

.04z 


1 
1 ft 


no 
Uo 


J 1 


^A 
.JO 


OO 

— 28 


07 

LI 


1 1 
1 1 


.Z 


n7 dt r 


c 

- J 


no 
llo 


.loooz 


KOI 
OOl 


1 7 


jy 


JO 


nn 
,UU 


OO 

— 28 


1A 
JO 


0/1 

z4. 


i 

. J 


1 1 "3 A 

1 1 3-A 




1 1 Q 

llo 


i o.o.n/1 
,looU4. 


i n^i 
lUol 


1 ft 


nn 
UU 




1. A 

. jO 


OO 

— 2y 


i 1 
1 1 


A/1 
U4. 


o 

-Z 


QC DT 


A 

■4 


llo 


i m no 

.lyiyz. 


ODD 


1 O 
1 ft 


n i 
UJ 


Uo 


nA 
.vo 


on 

— 2y 


1 
1 


Cl 

J J. 


.0 






11/1 
114 


1 OV1 

,iy ( iz. 


01 Q 
Olo 


1 
1 ft 


no 
Uz 


1 1 


^0 
.JZ 


OO 

— 2y 


1 o 

i y 


on 

ZU. 


A 
.0 


OA rt 


1 Q 

■iy 


11/1 
1 14 


1 O0./1 

■ iyo4o. 


111 
tit 


1 
1 ft 


no 
Uz 


JO 


70 

. /Z 


OO 

— 2y 


m 
UJ 


/i i 
4 1 . 


Q 

y 


11/1 A 




11/1 
114. 


i nn?n 

,iyy ( u 


O A Q 

o4o 


1 Q 

lo 


no 
Uz 


j4 


A Q 
.4ft 


on 

— 2y 


OA 

zo 


OO 

zy 


.4 


nc T3T 


0/1 

-Z4 


1 1 Q 

llo 


1 om a 
.loU14. 


oon 
ozU 


1 1 
1 / 


5 

jft 


5 

ZJ. 


on 
.zU 


on 

— 2y 


A 7 

4 / 


CO 

jy. 


.0 


n/i DT Z" 1 

y4-JrjLO- 


■3 


1 io 

llo 


1 Q 1 A 1 

-lol41. 


7Q 1 

. / ol 


1 1 
1 / 


Jft 


JO 


70 


— oU 


no 
Uz 


i n 

iy. 


.3 


n/i DT Z" 1 


-4 


1 io 

llo 


1 Q071 
-lOZ ( 1. 


700 
. ( OO 


1 / 


jft 


JZ. 


<A 
.JO 


— oU 


no 
Uz 


no 

Uo. 


-> 

,z 


lion 

1 lo-U 




11Q 

llo 


i 0/ino 
.lo4Uz. 


4yo 


1 7 
1 / 


cn 

jV 


1 1 
1 J 


OO 

.yz 


OO 

— 2y 


j j 


^7 

jj 


.4 


1 1 c 




11C 

llo 


1 Q7H7 

.io (\) ( 


i on7 

ioy / 


1 

1 ft 


nn 
UU 


nA 
UO 


OA 

,yo 


on 

— 2y 


jo 


nA 
UO 


n 

.U 


m dt r 


OA 

■zo 


11C 

llo 


i ni qo 
.lyioz. 


om 

oyi 


1 

1 ft 


n i 
UJ 


nn 

uv 


O ,1 

,o4 


on 

— 2y 


Z A 
JO 


1 o 

iy 


n 

.U 


no dt 


■33 


118 


.19184. 


939 


18 


01 


10 


.32 


—29 


48 


55. 


.4 


99-BLG- 


42 


121 


.22032. 


133 


18 


07 


46 


.80 


-30 


39 


41. 


.4 






158 


.27444. 


129 


18 


20 


10 


.32 


-25 


08 


25. 


.8 






162 


.25865 


442 


18 


16 


30 


.72 


-26 


26 


58. 


.2 






162 


.25868. 


405 


18 


16 


46 


.08 


-26 


11 


43. 


.1 


95-BLG- 


-8 


176 


.18826 


909 


18 


00 


15 


.36 


-27 


42 


47 


.9 






178 


.23531. 


931 


18 


11 


17 


.52 


-25 


59 


51 


.4 


178-A 




180 


.22240. 


202 


18 


07 


57 


.84 


-25 


24 


48. 


.2 






401 


.47991 


.1840 


17 


57 


07 


.44 


-28 


13 


44 


,0 


97-BLG- 


19 


401 


.47994 


1182 


17 


57 


07 


.92 


-28 


00 


27. 


.7 


98-BLG- 


■26 


401 


.48052. 


861 


17 


57 


23 


.76 


-28 


10 


24 


.6 


98-BLG- 


10 


401 


.48167 


1934 


17 


57 


57 


.36 


-28 


27 


09 


.7 






401 


.48229. 


760 


17 


58 


13 


.44 


-28 


20 


32. 


.3 


99-BLG- 


■25 


401 


.48408. 


.649 


17 


59 


08 


.88 


-28 


24 


54 


.7 


98-BLG- 


■14 


401 


.48469. 


.789 


17 


59 


22 


.32 


-28 


20 


21 


.8 






402 


,47678. 


1666 


17 


55 


22 


.56 


-29 


04 


12 


.7 


98-BLG- 


■39 



#16, BLG-12 

#11,BLG-11 BUL.SC35- 144974 
#7, BLG-13 1998-BUL-23 



#7, BLG-13 
#4, BLG-28 



1998-BUL-23 



#14, BLG-4 

1998-BUL-22 
#12, BLG-5 1999-BUL-07 



BUL.SC20-395103 
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TABLE 4 — Continued 



field.tile.seq RA (J2000) DEC (J2000) MACHO Alert EROS OGLE 



402.47737.1590 


17 55 41.52 


-29 09 00.7 






402.47742.3318 


17 55 42.96 


-28 49 17.4 






402.47796.1893 


17 56 10.80 


-29 10 44.4 


402-B 




402.47798.1259 


17 56 11.04 


-29 06 14.0 






402.47799.1736 


17 56 07.20 


-28 58 32.5 


98-BLG-37 




402.47856.561 


17 56 24.96 


-29 12 56.5 






402.47862.1576 


17 56 20.64 


-28 47 42.0 


97-BLG-41 




402.48158.1296 


17 57 47.52 


-29 03 54.4 


98-BLG-19 


1998-BUL-17 


402.48280.502 


17 58 24.96 


-28 57 46.4 






403.47491.770 


17 54 38.64 


-29 33 13.0 


99-BLG-7 


1999-BUL-02 


403.47550.807 


17 55 00.00 


-29 35 03.8 


99-BLG-54 


1999-BUL-41 


403.47610.576 


17 55 17.04 


-29 37 40.8 






403.47845.495 


17 56 22.80 


-29 55 16.3 


97-BLG-17 
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Table 5 

Events with light curve deviations 



field. tile. seq 


Comments 


101.21689.315 


possible parallax 


104.20119.6312/109.20119.1051 


parallax or xallarap 


104.20251.50 


parallax 


104.20382.803 


parallax 


104.20640.8423 




105.21813.2516 


possible red blend 


108.18951.593 


binary 


114.19970.843 


weak signal 


118.18141.731 


binary 


118.18797.1397 


asymmetry in blue peak 


158.27444.129 




180.22240.202 


post-peak blue points preferentially low 


401.48229.760 


binary 


401.48408.649 


binary 


402.47742.3318 




402.47799.1736 


coherent steep variation in blue 


402.47862.1576 


binary 


402.48280.502 


possible blue blend 


403.47610.576 


possible red blend 
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Table 6 

Event parameters from blend fits 
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^l.zoztU.zUJ X 1U 


( 1 nq 4- n 1 1\ i n — 5 
^l.Uo zt U.l / J X 1U 


(a nn 4- 1 en^ \^ i n — 6 
^y.uu zt l.DUj X 1U 


113 


.18804 


.1061 


104-1/1 

iy zt 14 


n 174- n ift 
U.l / zt U.10 


n 1 q 4- n i /i 

U.lo zt U.14 


n 1 a 4- n 1 e 

U.14 zt U. ID 


113 


.19192 


,365 


on 1 1 Q 
oU zt lo 


n Qft 4- n Q/i 

U.oO zt U.04 


1 K 4- o n 

i.d zt z.y 


o n 4- q 7 
z.u zt o. / 


114. 


.19712 


.813 


O/i n -U n 
z4.U zt z.U 


n oqq 4- n OQft 
U.zoo zt U.UoO 


n oq 4- n 1 ft 
u.yo zt u.io 


n 07 4- n 1 ft 
u.y / zt u. io 


114 


.19846 


.777 


1 OO -U 1 /l 

izz zt 14 


n OQ/i 4- n n^o 
U.zo4 zt U.UDz 


n ooft 4- n nr;i 
U.zzO zt U.UDl 


noii 4- n n/i q 

U.zll zt U.U40 


114. 


1 9970 

. -L y y t vj 


843 

.04:0 


1 ft Q _|_ 1 C 

lO.o zt 1.0 


n 1 01 4_n noc. 
U. lol zt U.UzD 


n 70 4- n i o 
u. / y zt u.iz 


n 7q 4- n 1 1 
U. / o zt U.ll 


1 1 8 


.18014 


390 


1 K 7 4- Q ft 

id. / zt 0.0 


1 ft _l_ 1 q 
1 .0 zt l.o 


C J_ 1 1 

D zt 1 1 


e 4- 1 1 

D zt 1 1 


1 1 8 


1 8971 


738 


OQ _|_ 17 

yo zt 1 < 


n oiq 4. n n^7 
U.zlo zt U.UD ( 


n qt; i n 1 1 
U.oO zt U.ll 


n qq 4_ n in 
U.oo zt U.1U 


1 1 8 


.18402 


495 


Oft 1 1 Q 

zO zt l.o 


n O/i ft 4- n no7 
U.z40 zt U.Uz ( 


n Qq i n 1 1 
U.oo zt U.ll 


n oc i n 1 1 
U.OD zt U.ll 


118. 


,18797 


.1397 


101 Q 4- O ^ 
IZl.o zt z.D 


n 1 07/1 4- n nnqo 
U.lz /4 zt u.uuoy 


1 nftn 4- n nq/i 
i .uoy zt U.Uo4 


1 c\a q 4- n nq/i 
l .U4o zt U.Uo4 


1 1 8 


,19182 


891 


10 /i q -U n 01 
lz.4o zt u.yi 


n 1 or; _i_ n nqn 
u. iyo zt u.uou 


i 01 4- n 1 o 
i.zi zt u.iy 


1 oo 4- n on 
l.zz zt u.zU 


118. 


,19184 


.939 


00 c _i_ re /I 
00. D zt D.4 


n Qfic 4- n n^7 

U.oOO zt U.Uo 1 


n ft7q 4- n nQ7 
u.o * o zt U.UO * 


n ft7e 4- n nQQ 

U.O ( D zt U.UOO 


1 91 
±z±. 


99n c t9 


■ loo 


1D.U zt y.o 


10 4-10 
1 .z zt 1 .z 


o q 4- k n 

z.o zt D.U 


O A 4- e O 
z.4 zt D.z 


158. 


,27444 


.129 


40. 1 zt o.U 


n Qfto 4- n nqft 
U.oOz zt U.UoO 


n 7i n 4- n nnq 
u. / iu zt u.uyo 


n ftQO 4- n non 
u.ooy zt u.uyu 


162. 


.25865 


.442 


1 i in 

bl ± 10 


0.75 ± 0.21 


1.01 ± 0.50 


0.95 ± 0.47 


162. 


.25868 


.405 


30.6 ±2.5 


0.386 ±0.057 


1.00 ±0.19 


1.04 ±0.20 


176. 


.18826 


.909 


3.68 ±0.38 


7.27 ±0.70 




1120 ±410 


178. 


.23531 


.931 


13.1 ±2.9 


0.54 ±0.22 


1.14±0.71 


1.11 ±0.70 


180. 


.22240 


.202 


454 ± 20 


0.183 ±0.010 


0.547 ±0.035 


0.551 ±0.035 


401. 


,47991 


.1840 


38 ±18 


0.77 ±0.62 


1.6 ±2.3 


1.6 ±2.4 


401. 


,47994 


,1182 


61 ± 12 


0.273 ±0.085 




0.97 ±0.36 


401. 


,48052 


.861 


83 ±15 


0.304 ±0.091 


0.74 ±0.28 


0.70 ±0.27 


401. 


,48167 


,1934 


48.5 ±8.2 


2.83 ±0.62 


37.1 ±24.0 


39.8 ±26.0 


401. 


,48469 


.789 


7.5 ±4.2 


< 0.26 


0.58 ±0.58 


0.54 ±0.55 


402. 


.47678 


,1666 


13.4 ±6.8 


0.45 ±0.43 


1.0 ± 1.4 


1.1 ± 1.5 


402. 


,47737 


.1590 


11.0±3.5 


0.111 ±0.055 


0.26 ±0.13 


0.29 ±0.15 


402. 


.47742 


.3318 


64.7 ±8.5 


0.065 ±0.020 


0.72 ±0.14 


0.61 ±0.12 
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TABLE 6 — Continued 



field. tile. seq 


i 






fb M 


402.47796.1893 


21.8 ± 1.4 


< 0.03 


0.801 ±0.081 


0.832 ±0.078 


402.47798.1259 


97800 ± 19500 


(6 ± 20219) x 10~ 7 


(2.59 ±0.52) x 10~ 4 


(2.31 ±0.47) x 10~ 4 


402.47799.1736 


80 ±44 


0.84 ±0.70 


2.5 ±4.0 


2.7 ±4.4 


402.47856.561 


52000 ± 16000 


(2.0 ± 1.1) x 10~ 5 


(3.8 ± 1.2) x 10" 5 


(5.0 ± 1.6) x 10" 5 


402.48158.1296 


20.7 ±4.9 


0.148 ±0.082 


0.56 ±0.26 


0.61 ±0.28 


402.48280.502 


119 ±28 


0.035 ±0.011 


0.064 ±0.018 


0.084 ±0.024 


403.47491.770 


45.3 ±5.3 


0.494 ±0.098 


1.19±0.35 


1.18 ±0.35 


403.47550.807 


10.9 ± 1.9 


0.33 ±0.12 


1.17±0.47 


1.13 ±0.45 


403.47610.576 


14.5 ± 1.2 


0.123 ±0.014 


0.492 ± 0.059 


0.608 ±0.073 


403.47845.495 


13.9 ±2.7 


0.62 ± 0.19 


1.23 ± 0.64 


1.20 ± 0.62 



NOTE. — Known binaries and two strong parallax events (104.20251.50 and 104.20382.803) are not included. 
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Table 7 

Known blend fractions for clump events with binary lenses 



Event 


frM 


fbM 


108.18951.593 (97-BLG-28) 


1.00 


0.90 


118.18141.731 (94-BLG-4) 


1.06 


1.04 


401.48408.649 (98-BLG-14) 


1.08 


1.10 


402.47862.1576 (97-BLG-41) 


0.95° 


0.92 a 




0.86 b 


0.84 b 



Note. — For the first three events, blend frac- 
tions are taken from Alcock et al. (2000c). The so- 
lutions for 97-BLG-41 are designated as follows: 
11 rotating binary solution obtained by Albrow et 
al. (2000). Johnson's I and V results are reported 
in f rM and fb M columns, respectively; 
b binary with a planet obtained by Bennett et al. 
(1999). 

For all cases, the statistical errors in blend fractions 
are < 0.05. 
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Table 8 
Data on the 94 bulge fields. 



FielH 
nciu 


I 


I, 


/v. /inn 


N 

1 v exposures 


fc ^UULleiy to j 


e (9nnHav«^ 


101 


3.73 


-3.02 


629.5 


804 


0.41 


0.61 


102 


3.77 


-4.11 


444 


421 


0.24 


0.46 


103 


4.31 


-4.62 


368 


331 


0.22 


0.37 


104 


3.11 


-3.01 


652 


1639 


0.40 


0.59 


105 


3.23 


-3.61 


539 


640 


0.39 


0.52 


106 


3.59 


-4.78 


346 


12 


< 0.003 


< 0.003 


107 


4.00 


-5.31 


269 


51 


0.007 


< 0.003 


108 


2.30 


-2.65 


790 


1031 


0.46 


0.74 


109 


2.45 


-3.20 


661 


761 


0.36 


0.59 


110 


2.81 


-4.48 


408 


650 


0.33 


0.55 


111 


2.99 


-5.14 


298 


305 


0.23 


0.37 


112 


3.40 


-5.53 


246.5 


43 


0.005 


< 0.003 


113 


1.63 


-2.78 


834 


1127 


0.49 


0.73 


114 


1.81 


-3.50 


617 


776 


0.37 


0.59 


115 


2.04 


-4.85 


325.5 


357 


0.24 


0.35 


116 


2.38 


-5.44 


268.5 


314 


0.21 


0.38 


117 


2.83 


-6.00 


198 


39 


0.006 


< 0.003 


118 


0.83 


-3.07 


741.5 


1053 


0.44 


0.74 


119 


1.07 


-3.83 


542.5 


1815 


0.45 


0.68 


120 


1.64 


-4.42 


394 


676 


0.35 


0.54 


121 


1.20 


-4.94 


344 


352 


0.22 


0.37 


122 


1.57 


-5.45 


266 


186 


0.15 


0.23 


123 


1.95 


-6.05 


210.5 


40 


0.006 


< 0.003 


124 


0.57 


-5.28 


265.5 


337 


0.22 


0.39 


125 


1.11 


-5.93 


189 


287 


0.21 


0.37 


126 


1.35 


-6.40 


182 


36 


< 0.003 


< 0.003 


127 


0.28 


-5.91 


192 


178 


0.13 


0.24 


128 


2.43 


-4.03 


516.5 


711 


0.35 


0.55 


129 


4.58 


-5.93 


234.5 


31 


0.006 


< 0.003 


130 


5.11 


-6.49 


169.5 


20 


0.02 


0.02 


131 


4.98 


-7.33 


130 


106 


0.12 


0.22 


132 


5.44 


-7.91 


102 


180 


0.17 


0.33 


133 


6.05 


-8.40 


70.5 


176 


0.15 


0.18 


134 


6.34 


-9.07 


46.5 


197 


0.09 


0.15 


135 


3.89 


-6.26 


179 


30 


0.004 


< 0.003 


136 


4.42 


-6.82 


137.5 


207 


0.17 


0.26 


137 


4.31 


-7.60 


101 


193 


0.15 


0.36 


138 


4.69 


-8.20 


95.5 


201 


0.16 


0.21 


139 


5.35 


-8.65 


81 


191 


0.11 


0.11 


140 


5.71 


-9.20 


64 


209 


0.13 


0.25 


141 


3.26 


-6.59 


159 


28 


0.03 


0.03 


142 


3.81 


-7.08 


126.5 


218 


0.15 


0.36 


143 


3.80 


-8.00 


97 


210 


0.17 


0.25 


144 


4.68 


-9.02 


50.5 


210 


0.10 


0.18 


145 


5.20 


-9.50 


49.5 


210 


0.12 


0.23 


146 


3.26 


-7.54 


106 


207 


0.19 


0.33 


147 


3.96 


-8.81 


66 


208 


0.19 


0.20 


148 


2.33 


-6.71 


161.5 


229 


0.18 


0.30 


149 


2.43 


-7.43 


112.5 


236 


0.19 


0.34 
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TABLE 8 — Continued 



Field 


/ 

I 


I, 




at , / 1 nn 

JVdumps/ 1UU 


AT 

t 'exposures 


e(50daysj 


e^zuuaays ) 


150 


2.96 


-8.01 


98.5 


229 


0.18 


0.28 


151 


3.17 


-8.89 


85 


194 


0.14 


0.23 


152 


1.76 


-7.07 


124 


235 


0.24 


0.32 


153 


2.11 


-7.87 


85.5 


231 


0.17 


0.35 


154 


2.16 


-8.51 


82 


247 


0.19 


0.20 


155 


1.01 


-7.44 


95.5 


247 


0.17 


0.36 


156 


1.34 


-8.12 


88.5 


249 


0.14 


0.33 


157 


0.08 


-7.76 


105.5 


258 


0.20 


0.37 


158 


7.08 


-4.44 


246.5 


260 


0.20 


0.29 


159 


6.35 


-4.40 


258.5 


464 


0.27 


0.47 


160 


6.84 


-5.04 


248 


208 


0.15 


0.26 


161 


5.56 


-4.01 


366 


467 


0.32 


0.45 


162 


5.64 


-4.62 


301 


382 


0.22 


0.38 


163 


5.98 


-5.22 


259 


251 


0.19 


0.36 


164 


6.51 


-5.90 


167.5 


20 


0.002 


< 0.003 


165 


7.01 


-6.38 


135 


14 


< 0.003 


< 0.003 


166 


7.10 


-7.07 


131.5 


110 


0.09 


0.14 


167 


4.88 


-4.21 


388.5 


364 


0.26 


0.40 


168 


5.01 


-4.92 


298 


210 


0.17 


0.26 


169 


5.40 


-5.63 


208 


24 


0.004 


< 0.003 


170 


5.81 


-6.20 


170 


17 


< 0.003 


< 0.003 


171 


6.42 


-6.65 


120.5 


125 


0.15 


0.19 


172 


6.82 


-7.61 


89 


132 


0.13 


0.18 


173 


6.92 


-8.41 


87 


146 


0.06 


0.09 


174 


5.71 


-6.92 


143.5 


144 


0.18 


0.22 


175 


6.10 


-7.55 


111.5 


97 


0.06 


0.10 


176 


2.93 


-2.30 


814.5 


423 


0.25 


0.43 


177 


6.75 


-3.82 


319 


416 


0.30 


0.43 


178 


5.24 


-3.42 


477 


376 


0.27 


0.41 


179 


4.92 


-2.83 


510 


349 


0.20 


0.39 


180 


5.93 


-2.69 


478 


343 


0.21 


0.38 


301 


18.77 


-2.05 




925 


0.31 


0.46 


302 


18.09 


-2.24 


- 


365 


0.28 


0.37 


303 


17.30 


-2.33 




369 


0.23 


0.41 


304 


9.07 


-2.70 




458 


0.30 


0.50 


305 


9.69 


-2.36 




411 


0.29 


0.46 


306 


8.46 


-3.03 




423 


0.28 


0.41 


307 


7.84 


-3.37 




435 


0.27 


0.43 


308 


10.02 


-2.98 




360 


0.26 


0.40 


309 


9.40 


-3.31 




347 


0.23 


0.40 


310 


8.79 


-3.64 




387 


0.26 


0.43 


311 


8.17 


-3.97 




396 


0.28 


0.40 


401 


2.02 


-1.93 


964.5 


429 


0.25 


0.45 


402 


1.27 


-2.09 


1153 


1256 


0.26 


0.45 


403 


0.55 


-2.32 


1085 


473 


0.26 


0.45 



NOTE. — Efficiencies (columns 6 and 7) are averaged over A, 
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Table 9 
Optical Depth by Field 



Field 


(l,b) 


A^ars/100 


Nclumps/100 


JV, events 


t/10- 6 


101 


(3.73,-3.02) 


6284 


629.5 


l 


1.62 ± 1.62 


102 


(3.77, -4.11) 


6545 


444 


l 


3 67 ± 3 67 


104 


(3.11,-3.01) 


5813 


652 


9 


8.76 ± 3.73 


105 


(3.23,-3.61) 


6375 


539 


1 


0.35 ± 0.35 


108* 


(2.30,-2.65) 


6498 


790 


6 


2.04 ± 0.92 


109* 


(2.45,-3.20) 


6926 


661 


2 


0.58 ± 0.41 


110 


(2.81,-4.48) 


6649 


408 


1 


0.46 ± 0.46 


113* 


(1.63,-2.78) 


6252 


834 


3 


0.55 ± 0.35 


114* 


(1.81,-3.50) 


6665 


617 


3 


1.19 ±0.74 


118* 


(0.83,-3.07) 


6347 


741.5 


7 


2.85 ± 1.35 


119* 


(1.07,-3.83) 


7454 


542.5 







121 


(1.20,-4.94) 


5855 


344 


1 


0.59 ± 0.59 


158 


(7.08,-4.44) 


4703 


246.5 


1 


2.96 ± 2.96 


162 


(5.64,-4.62) 


5914 


301 


2 


3.11 ±2.21 


176 


(2.93,-2.30) 


7741 


814.5 


1 


0.37 ±0.37 


178 


(5.24,-3.42) 


8186 


477 


1 


0.43 ± 0.43 


180 


(5.93,-2.69) 


6388 


478 


1 


4.62 ± 4.62 


401* 


(2.02,-1.92) 


6630 


964.5 


7 


5.13 ±2.16 


402* 


(1.27,-2.09) 


7098 


1153 


10 


3.95 ± 1.50 


403* 


(0.55,-2.32) 


6798 


1085 


4 


1.16 ±0.66 


CGR a 


(1.50,-2.68) 


60668 


7388.5 


42 


9 i 7 +0.47 
' -0.38 



NOTE. — Only fields with at least one clump giant event (and field 119) are shown. 
* Field is near the Galactic center and included in CGR. 
a Average towards the group of 9 fields designated as CGR. 
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Table 10 

one sigma upper limits on the optical depth for fields with no clump events 



field 


^_lim,ler 


Field 


^_lim,l(T 


Field 


^_lim,l(T 


Field 


^_lim,l(T 



103 


3 


26 


X 


10" 


6 


128 


1 


45 


X 


10 


-6 


144 


4 


62 


X 


10 


-5 


161 


2 


52 


X 


10 


-6 


106 


7 


90 


X 


10" 


6 


129 


9 


87 


X 


10 


-5 


145 


4 


10 


X 


10 


-5 


163 


5 


84 


X 


10 


-6 


107 


4 


22 


X 


10" 


5 


130 


1 


13 


X 


10 


-4 


146 


1 


42 


X 


10 


-5 


164 


3 


68 


X 


10 


-4 


111 


4 


■56 


X 


10" 


6 


131 


2 


01 


X 


10 


-5 


147 


3 


10 


X 


10 


-5 


165 


6 


81 


X 


10 


-4 


112 


7 


12 


X 


10" 


5 


132 


2 


04 


X 


10 


-5 


148 


9 


71 


X 


10 


-6 


166 


2 


52 


X 


10 


-5 


115 


3 


■58 


X 


10~ 


6 


133 


3 


25 


X 


10 


-5 


149 


1 


39 


X 


10 


-5 


167 


2 


83 


X 


10 


-6 


116 


5 


02 


X 


10~ 


6 


134 


5 


45 


X 


10 


-5 


150 


1 


59 


X 


10 


-5 


168 


6 


30 


X 


10 


-6 


117 


8 


37 


X 


io- 


5 


135 


9 


14 


X 


10 


-5 


151 


2 


67 


X 


10 


-5 


169 


1 


05 


X 


10 


-4 


119 


1 


10 


X 


io- 


6 


136 


1 


23 


X 


10 


-5 


152 


1 


18 


X 


10 


-5 


170 


4 


54 


X 


10 


-4 


120 


1 


91 


X 


io- 


6 


137 


1 


72 


X 


10 


-5 


153 


1 


83 


X 


10 


-5 


171 


2 


08 


X 


10 


-5 


122 


6 


90 


X 


io- 


6 


138 


1 


97 


X 


10 


-5 


154 


2 


14 


X 


10 


-5 


172 


2 


91 


X 


10 


-5 


123 


7 


70 


X 


io- 


5 


139 


2 


80 


X 


10 


-5 


155 


1 


65 


X 


10 


-5 


173 


4 


44 


X 


10 


-5 


124 


4 


94 


X 


io- 


6 


140 


3 


42 


X 


10 


-5 


156 


1 


85 


X 


10 


-5 


174 


1 


30 


X 


10 


-5 


125 


7 


95 


X 


io- 


6 


141 


9 


01 


X 


10 


-5 


157 


1 


39 


X 


10 


-5 


175 


3 


58 


X 


10 


-5 


126 


1 


35 


X 


io- 


4 


142 


1 


28 


X 


10 


-5 


159 


3 


57 


X 


10 


-6 


177 


3 


09 


X 


10 


-6 


127 


1 


02 


X 


10- 


5 


143 


1 


80 


X 


10 


-5 


160 


8 


12 


X 


10 


-6 


179 


2 


01 


X 


10 


-6 
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Table 1 1 

Events from the verification sample 





Event 




(5bl 






Event 








108 


18947 


3618 





04 


102 


22466 


140 





55 


1 AO 

108 


1 one 1 
18951 


1 OO 1 


a 
U 


31 


1 A /I 

104 


19992 


o c o 
8j8 


n 
U 


OA 


108 


18952 


941 





74 


104 


20119 


6312 


-1 


14 


108 


19334 


1583 





56 


104 


20259 


572 





02 


109 


20640 


360 





25 


104 


20515 


498 





66 


113 


19192 


365 





27 


104 


20645 


3129 


1 


08 


114 


19712 


813 


-0 


20 


105 


21813 


2516 





15 


118 


18014 


320 





40 


110 


22455 


842 





70 


118 


18402 


495 


-1 


35 


121 


22032 


133 





26 


118 


18797 


1397 


1 


28 


162 


25865 


442 


-0 


12 


118 


19182 


891 


1 


11 


162 


25868 


405 





21 


401 


47991 


1840 





27 


178 


23531 


931 





16 


401 


47994 


1182 


-0 


08 












401 


48052 


861 


-1 


13 












401 


48167 


1934 


1 


49 












401 


48469 


789 


-0 


83 












402 


47678 


1666 





07 












402 


47799 


1736 





40 












402 


48158 


1296 


-1 


39 












403 


47491 


770 





51 












403 


47550 


807 





28 












403 


47845 


495 





32 













NOTE. — Events in CGR are displayed in the left col- 
umn and events in the remaining fields in the right col- 
umn. 
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Table 12 

Optical Depth by Field from Verification Sample 



Field 


N vcr 
events 


r ver/ 10 -6 


AJ-orig 
events 


t/10~ 6 


108 


4 


2.25 ± 1.20 


6 


2.04 ±0.92 


109 


1 


0.59 ±0.59 


2 


0.58 ±0.41 


113 


1 


0.80 ±0.80 


3 


0.55 ±0.35 


114 


1 


1.07 ± 1.07 


3 


1.19 ±0.74 


118 


4 


2.77 ± 1.85 


7 


2.85 ± 1.34 


119 












401 


5 


5.41 ±2.67 


7 


5.13 ±2.16 


402 


3 


4.55 ±3.42 


10 


3.95 ± 1.50 


403 


3 


1.29 ±0.85 


4 


1.16 ±0.66 


101 









1.62 ± 1.62 


102 


1 


3.67 ±3.67 




3.67 ±3.67 


104 


5 


6.07 ±3.18 




8.76 ±3.73 


105 


1 


0.35 ±0.35 




0.35 ±0.35 


110 


1 


0.46 ±0.46 




0.46 ±0.46 


121 


1 


0.59 ±0.59 




0.59 ±0.59 


158 









2.96 ±2.96 


162 


2 


3.11 ±2.21 




3.11 ±2.21 


176 









0.37 ±0.37 


178 


1 


0.43 ±0.43 




0.43 ±0.43 


180 









4.62 ±4.62 
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